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Abstract 

The deformation equation of a spacelike submanifold with an arbitrary 
codimension is given by a general construction without using local frames. 
In the case of codimension- 1, this equation reduces to the evolution equa- 
tion of the extrinsic curvature of a spacelike hypersurface. In the more 
interesting case of codimension-2, after selecting a local null frame, this de- 
formation equation reduces to the well known (cross) focusing equations. 
We show how the thermodynamics of trapping horizons is related to these 
deformation equations in two different formalisms: with and without in- 
troducing quasilocal energy. In the formalism with the quasilocal energy, 
the Hawking mass in four dimension is generalized to higher dimension, 
and it is found that the deformation of this energy inside a marginal sur- 
face can be also decomposed into the contributions from matter fields and 
gravitational radiation as in the four dimension. In the formalism with- 
out the quasilocal energy, we generalize the definition of slowly evolving 
future outer trapping horizons proposed by Booth to past trapping hori- 
zons. The dynamics of the trapping horizons in FLRW universe is given 
as an example. Especially, the slowly evolving past trapping horizon in 
the FLRW universe has close relation to the scenario of slow-roll infiation. 
Up to the second order of the slowly evolving parameter in this generaliza- 
tion, the temperature (surface gravity) associated with the slowly evolving 
trapping horizon in the FLRW universe is essentially the same as the one 
defined by using the quasilocal energy. 
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1 Introduction 



Quantum mechanics together with general relativity (GR) predicts that black 
hole behaves like a black body, emitting thermal radiations, with a tempera- 
ture proportional to its surface gravity at the black hole horizon and with an 
entropy proportional to its horizon area [U [2] . The Hawking temperature and 
the horizon entropy together with the black hole mass obey the first law of 
black hole thermodynamics [3 . Since these seminal works in the 1970s, the 
relation among thermodynamics, quantum theory and spacetime geometry has 
been widely discussed, and recent status can be found in a nice review [4j. 

Most of the studies of the black hole thermodynamics have been focused 
on the event horizons of stationary black holes. For example, Kerr-Newmann 
solution family in GR. However, this kind of horizon strongly depends on the 
global structure of the spacetime. For example, to define an event horizon, we 
have to know the future null infinity of the spacetime. The so called apparent 
horizon defined by Hawking [5^ also depends on the slicing (or 3 + 1 decompo- 
sition) of the spacetime. To classify the two surfaces embedded in some slice 
(for example, trapped, untrapped and marginal trapped) of the spacetime, of 
course, one has to study the extrinsic properties of the two surfaces embedded 
in this slice. To define the apparent horizon (a hypersurface) , one has to repeat 
this classification in each slice of the spacetimes. Recent years, based on the 
Hawking's definition of the apparent horizon, people have given some defini- 
tions of the so called quasilocal horizons, see a review [6J. Further, from a more 
broad view, some local defined horizon has also been proposed by Jacobson et 
al, and profound connection between gravitation and thermodynamics has been 
revealed |7, 8_. 

The pioneer work on the quasilocal horizon is the trapping horizon defined 
by Hayward more than ten years ago [SJ [TUl [H] . Roughly speaking, this kind 
of horizon is a hypersurface foliated by marginal surfaces of the spacetime. 
However, here, the so called marginal is different from the one given by Hawking: 
The former is a two dimension surface embedded in the spacetime [121 , while the 
later is a kind of two dimension surface embedded in some slice of the spacetime. 
The trapping horizon can be null, spacelike or timelike according to different 
spacetime structures. At the end of 1998, Ashtekar et al proposed a new horizon 
which is called isolated horizon [131 HH HSl HSl [13 • This kind of horizon is null, 
and does not exchange matter and energy with it's environment. About two 
years later, dynamical horizon is also proposed by Ashtekar and Krishnan jl81 
[T^ . It's a spacelike hypersurface and dynamically evolves. Another interest 
object is the so called slowly evolving horizon developed by Booth et al |20[ 1211 
I22| . In this proposal, the marginal surface (and associated physical quantities) 
slowly evolves on the future trapping horizon (In fact, in this theory, one can 
also consider the event horizon [22)). The slowly evolving horizon describes 
some near equilibrium state of the thermodynamics of the spacetime. 

To study the dynamics of the quasilocal horizon, usually, there are two 
formalisms: One of them heavily depend on some quasilocal energy inside a 
given two dimension surface. By using Einstein equation and the quasilocal en- 
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ergy, one may directly gets some first law like equation. Actually, to study the 
dynamics of the trapping horizon in general spherically symmetric spacetime, 
Hayward has translated the Einstein equation into a simple first law like equa- 
tion by selecting the quasilocal energy to be Misner- Sharp energy [3^. Without 
the spherical symmetry, the problem becomes complicated. Based on the Hawk- 
ing mass (energy) [^7] and focusing and cross focusing equations, the dynamics 
of the trapping horizon is also studied by Hayward et al [28l [29l [30l [31] . In 
this formalism, the quasilocal energy plays a key role, while the (cross) focusing 
equations link the variation of the quasilocal energy, matter fields and some pos- 
sible gravitational radiation together. Another formalism is independent of the 
quasilocal energy. In some sense, the most important implement in this method 
is the (cross) focusing equations (or generalized version]^. From these focusing 
(cross focusing) equations, firstly, one can defines some energy flux which is nat- 
urally related to the variation of the area of the cross section of the quasilocal 
horizon. This is a Clausius like equation if we simply regard that the area of the 
cross section corresponds to the entropy associated with the horizon. Secondly, 
one can also study the variation of the angular momentum of the horizon. Fur- 
ther, by assuming the first law of thermodynamics is still valid on the quasilocal 
horizon, from the Clausiu like equation and the variation of the angular momen- 
tum, one finally gets some energy of the horizon. In this formalism, we need 
not introduce some quasilocal energy in advance. Contrarily, the energy of the 
horizon can be regarded as a byproduct of the theory [SI [THl [TIJ [5D1 [HI |25] . 
However, these two formalism are both useless for the isolated horizon. In fact, 
there is no dynamical version first law of thermodynamics (or dynamical version 
Clausius relation) associated with the isolated horizon. To study this object, 
one has to consider phase space method [6l [131 [3 [El [El [E] ■ 

So the (cross) focusing equations are very important to study the dynamics 
of the quasilocal horizons. Actually, to study the dynamics of local Rindler 
horizons, Jacobson et al also apply the so called Raychaudhuri equation (cor- 
responding to one of the focusing equations) to relate the variation of the area 
of the horizon and the matter flux [71 [S] . Therefore, these equations inevitably 
appear whenever we hope to study the dynamics of the horizon quasilocally or 
locally. On the other hand, the (cross) focusing equations can be get from the 
study of the deformation of the codimension-2 spacelike submanifold. In the 
light of the importance of these equations, in this paper, we study the defor- 
mation of a spacelike submanifold with an arbitrary codimension, and a local 
frame independent equation will be given. In the case of codimension- 1, this 
equation reduces to the evolution equation of the extrinsic curvature of a space- 
like hypersurface in an n-dimensional Einstein theory. In the more interesting 
case of codimension-2, after selecting a local null frame, this equation naturally 
reduces to the well known (cross) focusing equations. 

According to the two formalisms mentioned in previous paragraph, in this 
paper, we show how the thermodynamics of the trapping horizons are related 
to these deformation equations in two different approaches: with and without 

^The Einstein equation is always used in any formalism. 



3 



introducing the quasilocal energy. In the first formahsm, we generahze the 
Hawking mass (energy) in four dimension to higher dimension and study the 
deformation of this energy. We find the deformation of this energy inside a 
marginal surface (More precisely, the evolution of the energy inside the marginal 
surface on the horizon) can be also decomposed into the contributions from the 
matter fields and the gravitational radiation as in the case of the four dimension. 
When the marginal surfaces are closed Einstein manifolds, we also study the first 
law like equation of the trapping horizon. Roughly speaking, it has a similar 
form as the one with the spherical symmetry. However, generally, it's impossible 
to define a surface gravity which is a constant on the marginal surface. Further, 
the surface gravity also evolves on the trapping horizon even in the spherically 
symmetric case. This means the system is generally nonequilibrium (even far 
from some equilibrium point) if we regard the temperature is proportional to 
the surface gravity. To make the problem easy to understand, it's necessary 
to study some near equilibrium state at first. The slowly evolving horizon 
is just the object which describes the near equilibrium state of the horizon 
thermodynamics. However, most of the studies of the slowly evolving horizon 
have been focused on the future outer trapping horizons [20l [211 HI], and it's 
quite necessary to study a slowly evolving past trapping horizon if we hope to 
study the near equilibrium state of some spacetime which only contains the past 
trapping horizon (for example, our cosmology). Therefore, in the formalism 
without introducing the quasilocal energy, we generalize the definition of the 
slowly evolving future outer trapping horizon to the past trapping horizon. 

The dynamics of the trapping horizon in FLRW universe is studied as an 
example. We find the slowly evolving past trapping horizon in the FLRW uni- 
verse has close relation to the scenario of the slow-roll inflation. With the slowly 
evolving conditions, we also find: up to the second order of the corresponding 
slowly evolving parameter, the temperature (surface gravity) associated with 
the slowly evolving past trapping horizon in the FLRW universe is essentially 
the same of the temperature (surface gravity) defined by using the quasilocal 
energy. 

This paper is organized as follows: In Sec. 2, a preliminary of submanifold 
theory is given. In Sec. 3, we deduce the deformation equation of a spacelike 
submanifold with an arbitrary codimension, and the cases of codimension-1 
and codimension-2 are given as examples. In Sec. 4, By selecting some local 
frame, we show our equation reduces to the (cross) focusing equations (in the 
case of codimension-2). We also give the deformation equation of the iS'0(l, 1) 
connection associated with the local frame. The general definition of the trap- 
ping horizon is provided in Sec. 5. In Sec. 6, we generalize the Hawking mass of 
the four dimension to the higher dimension, and study the deformation of this 
quasilocal energy by using the general deformation equation we have obtained. 
The dynamics of the trapping horizon is discussed based on the deformation of 
this generalized Hawking energy. In Sec. 7, Firstly, we study the null trapping 
horizon which corresponds to some equilibrium state, and then we generalize 
the slowly evolving future outer trapping horizon to the past cases. In Sec. 8, 
the tapping horizon in the FLRW universe is investigated. We give the slowly 
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evolving conditions of the trapping horizon in the FLRW universe. The surface 
gravity associated with this kind of trapping horizon is also studied. The Sec. 9 
is conclusion and discussion. 

2 Theory of Submanifold and Geometry of Codi- 
mension -2 Surface 

To define the trapping horizon, one has to study the intrinsic and the extrinsic 
geometries of the codiniension-2 spacelike surface, S, embedded in the space- 
time {M,g), and then define the union of all the codimension-2 surfaces which 
satisfy some conditions (be marginal) to be the horizon of the spacetime. So 
it's important to study the geometry of the codimension-2 surface. 

To describe the geometry of the codimension-2 surface, in this section, we 
list some important formula in the submanifold theory. We use abstract index 
notation to make the formula have similar style as the theory of the hypersurface 
in GR [35] . More details can be found in the papers of Carter [331 IMl |3S] . The 
readers who are familiar with the submanifold theory can skip this subsection. 

For a spacelike codimension-2 surface, from the submanifold theory, one can 
always decompose the metric of the spacetime into 

gab = hab + Qab , (1) 

where qab is the induced metric of the surface S. The induced metric qab is 
Riemannian because that S is spacelike, while the transverse part, i.e., hab, is 
Lorentzian. By this decomposition, the corresponding projection operators are 
given by qj' and hj", and they satisfy 

QaQc'^Qa, h:K' = K\ h:qJ' = Q. (2) 

Since hab is Lorentzian, it's natural to introduce two future directed null vector 
fields £ and n, and express hab as 

hab = - naib = sijeiei . (3) 

where / and J take values {1, 2}, and = £, = n. The symbol £/j represents 
a constant matrix given by en = £22 = 0, £12 = £21 = —1. According to 
pointing to singularity or not, the vectors n and € are called inward (ingoing) 
and outward (outgoing) respectively. Obviously, there are some freedoms to 
choose £ and n. However, if we require ta^"" = — 1, the only remainder freedom 
is just the rescaling of the null vectors, i.e., t — > Ai?, n n/X with some positive 
regular function A. 

Certainly, one can also introduce an orthogonal frame such that hab can be 
expressed as 

hab = ~UaUb + VaVb , (4) 

where Ua and Va satisfy: u"'Ua = — 1 , v'^Va = 1 and M°Wa — 0. Now, in eq.Q, 
we can take — u, e'^ — v, en = —£22 = —1 and £12 = £21 = 0. Similar to the 
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null frame, there are also some freedoms (for example, the SO{l, 1) rotation of 
the frame) to take different Ua and Va- 

Assuming the covariant derivative of the spacetime (M^g) is given by V, 
then, second fundamental tensor K^^ is defined as |33) 

Kab"" = qaUb^dQe' ■ (5) 

This is an important extrinsic quantity of the surface S*, and it can be defined 
without introducing any local frame of the spacetime. It is easy to find that the 
second fundamental tensor satisfies 

Kab' - /^6a" , K^Kdb' = Kab%' = . (6) 

This tensor can be decomposed into a traceless part (C^f,*^) and a trace part 
(K-), i.e., 

Kab"" ^ ^-^QabK' + Cab" , (7) 

where K'^ = g°'^Kab'^ is called extrinsic curvature vector or mean curvature 
vector^ which is an important tensor in the submanifold theory. By using the 
null frame, one gets these extrinsic quantities along the directions of £ and n: 

K^b = -Kab'tc = qa'qb'^Vdd , K'^^b = ~Kab'n, = qa'qb'^V.Ud . (8) 

Similarly, the extrinsic vector is also decomposed as 

= = 9«''V„4 , e^"' - ~K'n^ - g-^'Vanb . (9) 

These two quantities are called the expansions along £ and n respectively. The 
traceless part is decomposed as 

all' = -Cab'lc = [qa'qb'' - ^^^labq""^ Vc4 , 

'y^ab = -Cab'n, = (^qa^qt,^ - V.rid . (10) 

These are just the usual shear tensors along the directions of t and n. For the 
orthogonal frame {u,v}, we can also get k'^^ = —K^j^Uc, Kj^J = —K^^Vc 
and the corresponding expansions and shear tensors. Actually, for an arbitrary 
normal vector X, we can define 

and the expansion and the shear tensor are respectively given by 

To study the intrinsic geometry of S*, it's necessary to introduce the corre- 
sponding connection or covariant derivative Da on S. For the tensor field which 
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is invariant under projection operator qj", for example, Tj', the corresponding 
covariant derivative is defined by 

DcT," = q^Ue\J^}T/ , (11) 

and more general cases are similar. With this definition of the covariant deriva- 
tive, for vectors ^ and r] which are tangent to S (invariant under the projection 
operator qj'), it's easy to find 

ry-Ve^" = 7J-D,e + KacVe ■ (12) 

This is just Gauss's formula. By using this covariant derivative Da, from the 
usual definition 

Rabcd^'^ = {DaDb - Dt,Da)^, , (13) 

one gets the intrinsic Riemann curvature tensor Rabcd- The relation between 
this intrinsic curvature of S and the curvature of the spacetime is encoded in 

Gauss equation: 

Rabcd = Kca^Kbde " K^b^Kade + Qaqb^Qc^Qd^ ^efgh , (14) 

where .'^abcd IS thc Riemann curvature of the spacetime. This equation can be 
easily found from the definitions of Da and Rabcd- 

Similar to the covariant Da for the intrinsic geometry of S, for a mixed 
tensor Tabc--- with tangent indices a, b and a normal index c, it's convenient to 
define a covariant derivative Da as follows: 

DeTabc- = /i/ • • • Qa'Qb" ■ ■ ■ Qj^jTghd ■ (15) 

Obviously, for tangent tensor (which is invariant under q^), this covariant 
derivative reduces to the derivative Da- For an arbitrary normal vector X 
and a tangent vector ^, by using above definition, it's easy to find 

(^"VcX^ = -K^cCX" + ^DcX'' . (16) 

This is just Weingarten's formula. So, for normal vectors. Da is just the usual 
normal covariant derivative. Based on this covariant derivative, by calculating 

^abcdX" = {DaDb - DbDa)Xe (17) 

for an arbitrary normal vector X, we get the corresponding curvature tensor 
^afccd) which has form 

^abcd = qa^q/hj^hj^^^efgh + KaedK,% - KbedKa% • (18) 

This is Ricci equation. Obviously, this curvature has property of Weyl tensor, 
in fact, after some rearrangement, it can be expressed as 

i^abcd = qal/K'K^'^efgh + CaedCb\ - CbedCa\ , (19) 
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where '^efgh is the Weyl tensor of the spacetime, while Cabc is the traceless part 
of the second fundamental tensor. 

Further, in our codimension-2 cases, from ea.(|15p. for an arbitrary normal 
vector Xa = a£a + (3na, it's easy to find 

DaXb - [Daa + Uaa)h + {Dal3 - UJal3)nb , (20) 

where uia is defined as 

LUa = -q^^ndVeH." . (21) 

This is just the normal covariant derivative given in some references (for exam- 
ple [22]). Sometime, the Wa is called the 5*0(1, 1) connection of the SO{\, 1) 
normal bundle (see, for example, [13]). Of course, this definition of the connec- 
tion on the normal bundle depends on the null frame (so it's gauge dependant). 
Similarly, if we consider the orthogonal frame for Xa — aua + Pva, it's 

easy to find 

baXb = {Daa + UJaP)ub + {DaP + t^aa)vb , (22) 

and now uja is defined as 

LOa = -qa^u'^VcVb ■ (23) 

It should be noted here: we have used the same notation Wa as in the case of 
the null frame, but their values are usually different from each other. 
Generally, the connection is defined to be 

Wabc = £ljeiDae( = UJa^bc , (24) 

where eab = "-a^b — (^a^b for the null frame {£,n}, and LOa is given in eq. (|2ip . 
While for the orthogonal frame {u,v}, eab = UaVb — VaUb and uJa can be found 
in eq. (The properties of the tensor eab is given in Appendix B.). It's easy 
find this uJabc satisfies standard relations 

i^abJ" = Dalb , i^abcn" = DaUb . (25) 

for the null frame, and 

i^abcU" = DaUb , IjJabcV" = DaVb ■ (26) 

for the orthogonal frame. For the null frame {i, n} or the orthogonal frame 
{u, v}, from ea. (ll7p and above two relations, it's easy to find that the curvature 
tensor ^abcd now can be put into 

^abcd = DaUJbcd — DbLOacd + '^ac^bed — ^bf^aed ■ (27) 

Considering the definition of Da, this curvature is just the one proposed by 
Carter [33 [Ml [35]: 

^abcd = (hfh/qJ^q^''VgLUhef + L^ac^bed^ - (a O 6) . (28) 
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Actually, the expressions ^7} and (^51) are valid in the case with an arbitrary 
codimension. In the special case of codimcnsion-2 in this paper, after a short 
calculation, we find 

^abcd = {DaUJb - DbUJa)ecd ■ (29) 

Here, flab = -Da^t — Dj^uja is the curvature associated with the 5*0(1, 1) connec- 
tion (pri) . 

Another important formula is Codazzi equation. This equation can be ob- 
tained from applying the covariant derivative Da on the second fundamental 
tensor. From equation 

bd{KabcX') = bd{Kabc)X' + KabcbdX' , (30) 

and the relation K^^^Xc ~ —q^qJ^eXf, it's not hard to find: 

DaKbcd - DbKacd = -qa<lb^qc''hf^efhg • (31) 

This is just the Codazzi equation. For an arbitrary normal vector Y, it gives 



+ ~ Ka\bbY'' = q^q^^^L 



(32) 



In the case of codimension- 1, this equation is just the so called momentum 
constraint equation in Hamiltonian formalism in GR if we select = tt°, 
where m° is the unit normal vector of some spacelike hypersurface. In the case 
of codimension-2, by using eq. (j25p . immediately, we get 

(j^) iDa - UJa) 0^'^ - [Db - UJb) Ct'^T = < ebcd (33) 

and 

+ """^ ~ + "^"^ " ^'t^n'^S^.^bcd ■ (34) 

The right hands of above equations can also be transformed into the form com- 
posed by the Weyl tensor 'loabcd and Einstein tensor '^ab'- 

q^q^^VSi^bcd = qa'q'^Y^^bcd " ( ^) q^Y^'^eb ■ (35) 
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We can get similar equations 
'n — 3 



n-2j (^"^^"' ~ "^"^^"0 " l^"'^^"^'' - ^bCT^a") = q^q^'W'^^.bcd , (36) 

^) (dJ^-'' - - [ObCT^T - ^b^^a^') - qa'q'^V^ebcd (37) 

by using eq. ((26|) if the orthogonal frame is considered. Eqs. dTi]) . (|18l) and (|3ip 
are important relations in the submanifold theory. They are valid in the case 
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with an arbitrary codimension. In the four dimension, some of the formula 
for the codimension-2 surfaces we have hsted here can also be found in [22] . 
It should be noted here: Some definition might be not well defined if we only 
consider a single surface. For example, Va is not well defined if the support 
of the tensor field is confined to the single surface. Actually, in this case, only 
la^b is well defined on the surface [31]. However, in this paper, we do not 
consider this possibility. In our setting, we can always imagine that 5 is a leaf 
of the foliation of some neighborhood of S. By this consideration, Va is always 
well defined. 

3 Deformation of Submanifold 

3.1 Cases with Arbitrary Codimensions 

Now, let's consider the deformation of some spacelike surface S embedded in the 
spacetime with an arbitrary codimension. Assume X"^ is a normal vector, then, 
generally speaking, the Lie derivative of the projection operator qj' along X" is 
not vanished. However, we can consider a foliation of some neighborhood of S 
(Certainly, S" is a leaf of this foliation), and require that X has some relation to 
the structure of this foliation such that the projection operator is Lie dragging 
along X, i.e.. 



where Cx is the Lie derivative along X. This relation just means: any tangent 
tensor (invariant under qJ') preserves to be a tangent tensor under the Lie 
derivative along X. For example, for any tangent vector we have £x£,°' = 
^x[S}'%°') = Ib'C.xS!'- So the resulting vector is still a tangent vector. The 
relation (|38p also means that 



This equation can be easily found from the relation qj^ + hj^ = Sj'. 

Of course, now X is constrained by above two equations (only one of them 
is independent). The situation is very similar to the 3 + 1 decomposition in GR. 
In that case, for some foliation parameter r, the evolution vector X is required 
to satisfy Cxt = 1, and the three dimension projection operator is also Lie 
dragging along the evolution vector X. Here, we have generalized this relation 
to the case with an arbitrary codimension. For the case of codimension- 1, more 
details can be found in reference [36] . 

Contrarily, for a given normal vector X, we can find an appropriate foli- 
ation of some neighborhood of the given codimension-2 surface such that the 
corresponding projection operator is Lie dragging along X. 

The relation (j38p automatically implies following two equations: 




(38) 




(39) 




(40) 



X'^hbcVdK' = hMqa'VcX'' = DaXb . 



(41) 
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So, for another normal vector Y (which is arbitrary), we arrive at 

q.^X'^VdY, = -YcDaX' . (42) 

This equation just provides the expression of the tangent part of the vector 
X'^WdYc- Generahy, it may has a normal part. The second equation in eg. ipH)) 
gives result 

qJ'X'^V.q,^ - D'Xa - -h^q^'V.X" (43) 

where we have used eq. (j4ip . So, by considering the definition of Da in eq. (jl5p . 
we have 

h^Ua'i'^cXd + VdX,) ^ h^Ua't^Xgcd = , (44) 

This relation just means: after the deformation, the spacetime metric still pre- 
serves the orthogonal property of the projection operators. Furthermore, for 
the (arbitrary) normal vector Y , ea. (|43l) also implies relation 

[X,r]^qe6 = 0, (45) 

This just tells us that the Lie bracket of X and Y is also a normal vector. 
Actually, if this condition is regarded as a primary assumption, then, one can 
easily get the equation (|42|) . This logic has been used in the paper by Booth [22] . 
Now, it's easy to find that equations 

q.Ua'Cxq," = h%'CxqJ' - qt'K^CxqJ" = K'^K^CxqJ' - (46) 

are trivially satisfied. So there are no further useful constraint conditions pro- 
vided by eg. ([55]). 

By using the property that projection operator is Lie dragging along X, we 
can get a simple expression for the second fundamental tensor: 

Cxqab = qa^Qb^'^xqcd - -2Ka,^X, = 2A'(f ) . (47) 

Similarly, one finds the expansion along X can be expressed as 

Cxeg = e^^'h, , (48) 

where is the area element of the {n — 2)-dimension submanifold 5. Actually, 
once some quantity is invariant under the projection operator, we can use (|38p 
to simplify the Lie derivative of the quantity as above two examples. 

Let's consider the deformation of the second fundamental tensor K^^^^ (pro- 
jecting along the direction of Y). From above discussion, we have 

^xK^^ ^q^q^LxK^^ g„ 9f, A VgA^^ + K'^ a Kcb +^ b ^ca ■ 

(49) 

After a little bit complicated calculation (especially for the first term in right 
hand of above equation), we gets 

^xKZ^ - q:q,''X'^Yf^%,,as + K^^^'Ki"^^ - Y^DaDbX, 

+Kacb (>d^"A'^) - Kabc {X'^VdY') . (50) 
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The details to get above equation are given in Appendix A. Considering 
(£x(?"^)<Zc6 = -2(?'^^i^if ' and q:CxKb = Q, 

we have 

= {i^xKZ^)g''' + KZ^Cx {q'^' + h'^') 

= {CxKZ^9^'-2kZ'k^^'^'. (51) 

After substituting this resuft into eg . ([SO)) , we get the deformation equation of 
the expansion 0^^^ (The normal vector Y is arbitrary.), 

Cxe^''^ = q'^X'^Yf^ecdf - ^(^''^''if^f - Y-DaD-X, - {X^VdY^) . (52) 

The resuhs (ISUl) and (|5^ are vahd in the case with an arbitrary codimen- 
sion. In fact, in above discussions, we have not impose any requirement on the 
dimension of the part associated with hab- 

In above discussions, we only consider the deformation of the submanifold 
along a normal vector X°-. For a tangent vector, for example, 0°, the Lie 
derivative of 6'^^^ along 0" is constrained by the Codazzi equations (|3T|) and 

HMD: 

(^),,,.v,.,.^.,,.._(li^),.,,,^„,.. 

+ rC^^b^Y" + qf^^^'Y'^^.fgh ■ (53) 

Eqs. ([5n|) and ([5^ are main results of this subsection. To understand them, in 
following two subsections, we will apply them to the cases of codimension-1 and 
codimension-2 respectively. 

3.2 Codimension-1 Cases 

In the case of codimension-1, we can set hab — ^UaUh, where it° is an unit 
timelike normal vector of the hypersurface (the observer associated with the 
vector is just the so called Euler observer). So the extrinsic curvature is 
simply given by Kabc — KatUc- In this case, X is just the evolution vector 
Xa — Nua with lapse function N. We can select Ya — Ua such that 6**^^^ is 
given by 6^^'^ = K = —K°'Ua- Remembering the definition of Da in eq.(fT5|). we 
have 

Y'DaD'^X, = -DaD'^N , K^y)-bK^^f) ^ NK'^'^Kab , 
qSax-Y'^^.fgh = -NMabu'^u'' , K, [X'-V dY") = . (54) 
So eQ. (|52p is transformed into a very familiar form 

- ^LxK = .^,abU''v' + K'^'Kab - ^D-DaN . (55) 
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This is nothing but the evolution equation of K in the 3 + 1 formuhsm of GR 
(Actually, now, it's (n — 1) + 1 decomposition of the n-dimensional Einstein 
theory). It should be noted here: Usually, in GR, people use the definition 
of the extrinsic curvature as Kab = —qa<lb''^cUd- So it has a sign difference 
from our definition. If we use this definition, the minus in front of CxK/N will 
disappear. Similarly, we can get the evolution equation of Kab from ea.([50| and 
the Gauss equation (fT4|) : 

- ^CxKab = -qaqb'^^cd + Rab + KKab - 2KacK,' - ^DaDbN . (56) 

From this simple example, it's very clear how that the vector X is adapted 
with the foliation of the neighborhood of the hypersurface such that is Lie 
dragging along it. Of course, for a different lapse function N (or a different 
evolution vector X), one has to consider a different foliation structure of the 
neighborhood of the hypersurface [36] . 

3.3 Codimension-2 Cases 

Let's consider the more interesting case of codimension-2. Noted that in this 
case we have 

qfaX''Y''^%,fgh = -X'Y'^^.h - hf^X'Y'^^.fgh 

- -X'Y'^^^.h - hf^h'^^efgh {XbY") . (57) 

However, from the Gauss equation p^ . we get 

^ - 2Mabh'"' + h'''h'"'^abcd = R + KabcK"'' " K'K, . (58) 

After substituting eas. ([57|) and ([55]) into eg. ([5^ . we find 

CxO^^'> = -'Sab {X^Y^ - h''^ {XeY")] - K^^^'^^'Klf'^ 

+ i (i? + KabcK''''' - K^Ka) ■ (XeF^) 

-Y'^DaD'^Xa ~ {X'^VdY') , (59) 

where ^^a& is the Einstein tensor of the spacetime. We can rearrange the terms 
in above equation, and transform it into another more covariant form 

+ i (i? - KabcK'^''' - K,K') ■ {X.Y'') 

-Y^D^b^X^ - [X^V^Y") . (60) 

This is the main result of this subsection. Ea. (|60)) is very important to study 
the thermodynamics of the horizons. Of course, here, this deformation equation 
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is valid for an arbitrary codimension-2 surface (not only for the case which has 
closed relation to the horizons). 

During above reductions, the local frames = {£, n} or = {u, v} are not 
necessary. However, to find the explicit expression for the term Y'^DcD'^Xe, it's 
better to introduce some local frame. After some calculation (see Appendix C), 
we find 

Y'^D.D^X, = eij (Y'D'D.X^) + 2eij (ij'eMY'e'^D.X^) 

+ D,Lu'={eMY''X'')+u:,u:^X,Y'), (61) 

where Y^ is the component of along the vector e^, i.e., Y^ = e^. Da is 
the covariant derivative on the codimension-2 surface. For the null frame {£, n}, 
the matrix e/j and tensor eat are given in eqs.(|3]) and (j24[) respectively, and the 
5*0(1, 1) connection uJa can be found in eq. (PT|) . While for the orthogonal frame 
{u,v}, Eij and eab are given just bellow eqs.Q and (|24t respectively, and the 
connection uja can be found in eq. (1231) . 

From above results, it is easy to find that the operator Cx satisfies Cx-\-Cz = 
Cx+z- However, due to the existence of the term Y^DcD'^Xe, it has no property 
as the usual Lie derivative. For example, the action of the operator on the scalar 
0{Y) generally has property 

qfx)0^^^ ¥^ fCxO'^'^K (62) 

where / is a function on the spacetime. This behavior comes from the fact that 
we require that the projection operator is Lie dragging along X. Of course, in 
the case that / is a constant on the codimension-2 surface, from eq. (|6ip . one 
finds the operator £x reduces to the usual Lie derivative. We will find this 
point in the spherically symmetric cases, for example, the FLRW universe. For 
Y vector, we have 

Cxe'^f'-^^fCxe^^^ + e^'-^Cxf, (63) 

This is because that the derivative of the Y variable is first order in eq. (|60| . To 
make the problem easy to understand, in next section, we give the explicit form 
of the deformation equation in some local frames. 

The property (p^ of the operator Cx is very similar to the operator "i5" 
defined by Andersson, Mars and Simon [24j[25]. Here, to compare this operator 
6 with our approach, we give (the notations are a little bit different) the original 
definition of S [5S]: Assume G / C i? and X be an arbitray normal vector 
at S. Let $x ■ S x I ^ A4 he a differential map such that for each fixed 
T & I, ^x{-, t) is an immersion and for fixed p & S, ^x{p, t) is a curve starting 
at p with tangent vector X{p). This way, a family of codimension-2 surfaces 
Sr = ^x{S,t) is defined. Assume Y^ be a nowhere zero normal vector of Sr 
which is differentiable with respect to r, let O'^^^^ be the expansion of Sr along 
Ft-, then, the variation of 6'^^-' is defined as Sxd'^^^ = 9t^^'-^^-' |r=o- This is a kind 
of geometric variation and we can also assume that X has a part which is 
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tangent to S. By using this definition, and assuming X = + A£ — Bn, one 
gets [25]: 

+ 1 ( - 2watj" + 2A.CJ" + R + KaK'' - 2^ahrn'') , (64) 

where we have chosen Y — £. uja is defined in eq. ((2T|) . As = D'^Da is the 
Laplacian on {S, Qab), and a = — na9r^?|T=o- If we change notations as: X ^ q, 
£-^l,n->- k/2, uj-^ s, e^^'> 9, KaK"" -H^, A b and B ^ u, then, §^ 
is just the original result in the Lemma 3.1 in [25j . 

Our approach is different from the definition of S: For an arbitrary normal 
vector X of S, we assume that there is a foliation of some neighborhood of S 
such that qj' is Lie dragging along X . So our operator Lx is just the usual 
Lie derivative constrained by By this definition, we can get the value of 

CxO^^^ on S for each X. However, it will be found soon in the next section that 
the result given by 5x is very similar to our result by Cx (see the discussion 
bellow eg . (p7|) ) . In our approach, the action of is encoded in eq. ([55)) . To 
get the result of 5x0^^\ we need not know the details of the map $x- A 
different X may correspond a different map <i>. Similarly, to get CxO^'^\ we 
also need not know the details of the foliation structure of some neighborhood 
of S, and a different X corresponds to a different foliation of a (maybe different) 
neighborhood of S. 



4 Deformation Equations with Local Frames 

4.1 Expressions in Null Frame 

4.1.1 Focusing and Cross Focusing Equations 



Eqs. ([59|) and ([60| in last section are independent of any local frame. In this 
subsection, by choosing the null frame {£,n}, we give two important examples. 
The components of the extrinsic curvature along these two null vectors can be 
found in eq.(l8]), and the corresponding expansions and shear tensors are defined 
in eqs.® and (fTOj) respectively. 

Firstly, by setting Ya = £a and Xa = A£a — Bua, we have X^Y"^ = B. So 
eqs.(|5Ul) and (|5Tt give result 



£jf 6l(^) = Kjf6l(^) - D^D'^B + 2lo''DcB - B 



-A 



ah 



1 



(65) 



Here, we have introduced an important quantity — "surface gravity" 



Kx = -n'^X'^Ve 



(66) 
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and then the last term in cq. lpO)) becomes 



(67) 



Since we have not discuss any horizon until now, so, in eq. (|66p . nx associated 
with the vector X is not the real surface gravity of the black hole horizon. It 
should be noted here: if we identify "kx" with "a" in dMl), then (5x6*'^) in ^ 
is the same of our result ([S5|) . 

Similarly, by setting Yq = Ua and Xa — Aia — Bua, we get 



B 



1 



(68) 



In the four dimension, eqs. (|65|) and ([SS)) are just the equations given by Booth [22] 
In the case where A = 1, B = {or A = 0, B = -1), eq.^ ( or eq.(|68l)) 
is just the so called focusing equation. In the case where A = 0, i? = — 1 (or 
A = 1,5 = 0), ea.(|65| (or eq. (|68t ) gives the cross focusing equation (see eqs. dMl) 
in next section) [28 l [29 l |30] . 



4.1.2 Y is dual to X 

The second interest example is the case where X^Y^ — 0. Under this require- 
ment, ea. (l60l) becomes very simple: 

CxO^'''' = - i'^ab + K,daK'\) X'^Y'' - Y'D.D'X, - K, (X^VeF=) . (69) 

For simplicity, we select X and Y as 

Xa = Ala - BUa , Y^ = Al^ + BUa , (70) 

where A and B are functions on the spacetime. It should be noted here: above 
expressions are not the most general forms of X and Y (which satisfy the relation 
XgY'^ ~ 0). Actually, here, X and Y are dual to each other, i.e., Ya — eabX^ . 
The tensor eab has been defined in eq. ([M)) . After a short calculation, ea. (|5T]) 
reduces to 

Y'^DcD'^Xe = DeiAD^B - BD^A - 2ABlu'') . (71) 
It's also easy to find 

Kc (X^VeY^) = -Kx0^^^ - e'^^^CxA - e^'^^CxB . (72) 

So the deformation equation of the expansion 9'^^^ becomes 

£;,0(^) = nxe^^^ - '^abX-Y' ~ aif - 

-DeiAD^B - BD^A - 2ABuj'') + O^'^^CxA + e^""^ CxB , (73) 
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where kx is defined in ea. (|5Sl) . In four dimension, by setting A — 1, this 
equation is just the one given in [22l [38l |39] . Consider ea.(|63|. we can translate 
above equation into 



(X 



n — 2 



+D^{AD''B - BD^A - 2ABuj<') + ACxO'^^'^ + BCxO^"^ , (74) 

For n = 4, this result can also be found in [331 [351 [32] ■ This is a very important 
equation to study the thermodynamics of the horizon. If the codimension-2 
surface is compact without boundary, then, it's easy to find 

, \aCx9'^'^^ + BCx9'^"'^] , (75) 

where is the area element of the codimension-2 surface. Here we have used 
the relation 0^^hq = Cxeq in eq. pS]) . 



4.1.3 Damour-Navier-Stokes like Equation 

The reason to select the null frame {£, n} is not simply to make the formula easy 
to understand. In fact, without the local null frame, some important physical 
quantities can not be defined. For example, Ua and kx, which have close relation 
to the angular momentum and the surface gravity of the horizon. So another 
important equation is the deformation equation of tOa- From the definition of 
uJa in ea. ([3T|) . by using the relations (1351) and ((32), it's easy to find 

jCx^a = Kj'^Dtie^'^Xd) + Da^x - ^Qa' X'^e'^^dbce ■ (76) 

We can also express the right hand of above equation by the Weyl tensor "^abcd 
and the Einstein tensor "i^ab- 

qa'X'^e^-^dbce = qaX^e^-^dbce " ^^q^e^'^X^^c ■ (77) 
From the generalized Codazzi equation (1331) and eg. ([35)1 . we have 

- -qJ'h''Y''%,,a ~ (^) 'iaY'^bc , (78) 
where Ya is the dual vector of X^,, i.e., Ya = e-abX'^ . It's easy to find that 
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always holds. Thus, by combining eas. ([7^ . (177)) and ([75)) . we have 



n - 2 



- Z^e^lT^' + K.baY^ + -Z„V^^^fc, . (79) 



The term K^DaY'^ can be calculated from eas. (|70)) and (PH)). and above equation 
becomes 



^n- 2^ 

q,'Y'%c-0^''>DaA-9^"'>DaB. 



(80) 



In the case where X is self-dual or anti-self-dual, i.e., X = ±Y, by consid- 
ering the Einstein equation, this equation is a kind of D amour- Navier- Stokes 
equation [SJ. In this case, X or F can be identified to be the evolution vec- 
tor of the event horizon of the spacetime (Actually, when X = ztF, eq. (|80)) 
can be explained to be the Damour-Navier-Stokes equation on an arbitrary null 
hypersurface of the spacetime). If eq. (|80)) is applied to the trapping horizons 
instead of the event horizon, one gets generalized Damour-Navier-Stokes equa- 
tions [SHI SO]. In this case, X is identified to be the evolution vector of the 
trapping horizon, and Y is just the normal vector of the trapping horizon. Cer- 
tainly, now, X need not to be self-dual or anti-self-dual, i.e., X need not to be 
null. More details on the physical meaning of the terms in ea. ()80p can be found 
in the papers by Gourgoulhon et.al. [551 HO) . 
Let (jf" be a tangent vector, then, we have 



rt — 3 



DaKx 



n~2 



qJ'Y-%c~0<^'^DaA~d^''^DaB 



(81) 



Assume that ^" also satisfies Lxtf)"" = and Da(p'^ — 0, then we get 
Cx [egir^a) = 



(82) 



Here, we have assumed that the codimension-2 surface is compact without 
boundary. It's obvious this equation does not depend on the dimension of 
the spacetime. The left hand of above equation has closed relation to the an- 
gular momentum of the spacetime |41| . Actually, the angular momentum can 
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be defined as = / eq((/)°Wa). From above equation, it's not hard to find the 
balance equation of the angular momentum |38j. 

Eqs. (j73p and (j80|) have important applications in the theory of membrane 
paradigm 021 HS] • Recent progress on this topic can be found in |38[ HSl |35] 
and references therein. 



4.2 Expressions in Orthogonal Frame 

Besides the null frame {£,n}, sometime, it's useful to express the deformation 
equations in the orthogonal frame {u,v} . For example, in the studies of the 
dynamical horizon, it's convenient to choose some orthogonal frame. Assume 
Xa — Aua + -Bwa, and Ya = Ua, then we have 

K,(X'V,Y'')^~e'^^'^Kx. (83) 

Similar to the "surface gravity" defined in eq. lpS)) . here, we have defined 

Kx = VcX^VeU' . (84) 

From eq. lpTj) . it's also easy to find 

VDcD^Xe = -D^DrA - 2w'=DcB - BD^w^ - Auj^uj'' . (85) 
Substituting above results into ea. (|60| . we get 

LxO^""^ = e^'^^x + D^DcA + 2uj''DcB + BD^uj" 



By similar calculation, for Ya = Va, we have 

CxO'^"'' = 6'(")kx - D^DcB - 2uj''DcA - AD^uj" 



(86) 



-A 



B 



ab 



UJcl^ + ^abU U 

{u)n(u) _ n{v}n{v) 



(87) 



Let Xa — Aua + Bva and Ya ~ Bua + Ava, then, Y is the dual vector of X, and 
they automatically satisfy X'^Ye — 0. It's easy to find 



Kx 



6l(-^) - Dc [AD^B - BD'A + lo" {A^ - B^)] 



(X) (Y)ab _ 



^'^abX^Y'^a\t>a 
+e^'''^CxB + e^''^CxA 
Considering the relation (|63l) . we have 



(88) 



CxO'^^'^ = ACxO^"^ + BCxO^""^ + e^'''>CxA + O^^'^CxB , 
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then, when the codimension-2 surface is compact without boundary, we get 

+ACxO^^^ + BCx 



n - 2 



(89) 



We can also get the deformation equation of uJa, which is given by 



Lxcoa = DaKx + (^) DaO^""^ " D^o^l^^ + K^DaY^ + q,^Y'%c ■ (90) 

However, it should be noted here: now, Ua and kx are given in eas. (|23l) and 
([5^ respectively. Although the definitions of Wq and nx are dependent of the 
selection of the local frames, the deformation equation of cJq has similar form in 
different frames if we give an appropriate definition of the "surface gravity" kx- 
Actually, from above equation or ea.(|79p. the combination {Cx^a — DaUx) is 
frame independent. 

From the assumption Xa = Aua + Bva, Ya = Bua + Ava, and eg. d^^ . we 
can easily get the explicit expression of the term KcDaY'^: 

K.DaY' = ^DaB9<^"^ - DaAd^^^ ~ e^'^^ua ■ (91) 

By this, the deformation equation (PH)) becomes 



+ q,!'Y'%c-0^''^DaB-d^^^DaA. (92) 

Similar to get ea. ((82|) . it's not hard to find: for Da^" = and Cx4>'^ = 0, we 
have 

Cx J irt^a) = J e,i^^{D-<j>'' + D'r) + ^Y^'^ab 

+A4>''Dae^'''^ + Bc/j^DaO^''^ I . (93) 

Certainly, to get above equation, compactness of the codimension-2 surface is 
also required, and this equation has closed relation to the balance equation of 
the angular momentum of the spacetime. 

At the end of this section, we give some discussions: 

(i). The so called (cross) focusing equations are just the special examples 
of eg. (1601) . Principally, for a pair of arbitrary normal vectors X and Y, we 
can get the value of LxO''^^ on some given codimension-2 surface according to 
eq. (|60)) . Of course, in this procedure, we have to consider a foliation of some 
neighborhood of the codimension-2 surface such that the projection operator is 
Lie dragging along X. 
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(ii) . By selecting different frames, the deformation equations of uJaS have 
similar forms as given in eqs. d79l) and (|90)) . Further, the combination CxOJa — 
DaKx is frame independent. To get this conclusion, we have to define an appro- 
priate "surface gravity" kx to match the definition of the SO{l, 1) connection 
oJa- When X is self-dual or anti-self-dual, i.e., Y = ±X, we get the Damour- 
Navier-Stokes equation. 

(iii) . The situation where Y is the dual of X is very simple and special. The 
compactness of the codimension-2 surface is important. Without this assump- 
tion, we can not get the simple expression for the integrals of the deformation 
equations, i.e. eas. ([75)) and (15^ (or eas. ([Ml) and ([M)) in the orthogonal frame). 

(iv) . Until now, we haven't discussed some special hypersurface of the space- 
time. In fact, we have focused on the codimension-2 surfaces. Of course the 
horizon of the spacetime is kind of hypersurface. In next section, we will find 
how to study this kind of hypersurface by our knowledge on the geometry of 
the codimension-2 surface. 

5 Trapping Horizon 

In black hole theory, one of the most important objects is the so called event 
horizon which is the boundary of the causal past of the future null infinity [5ll32]. 
So, to describe the event horizon, one has to know some global information of 
the spacetime, for example, the future null infinity. This kind of horizon is a 
null hypersurface of the spacetime. However, the so called trapping horizon is 
very different from the event horizon [21 [TOl E] ■ The definition of the trapping 
horizon is quasilocal, and it does not depend on the asymptotical behavior of 
the spacetime. This kind of horizon can be null, spacelike or timelike according 
to different spacetime structures which are involved. 

The codimension-2 spacelike surface with = is called marginal 

surface. The surface with > is called trapped, and 0(^)6'^"-' < is 

called untrapped. 

By the definitions of the trapped surface and the untrapped surface, we can 
define two regions in the spacetime: A trapped (untrapped) region is the union 
of all trapped (untrapped) surfaces. 

We can give similar definitions by using the extrinsic curvature vector if" 
from the relation K'^Kc — — 20^^^0^"). Sometime, the formulism without the null 
frame is enough. However, to give a detailed study of the marginal surfaces, 
for example, to give a classification of the marginal surfaces, it's inevitable to 
introduce the null frame or some similar structure. 

A marginal surface is called future if O'^^^ — 0, 6''^"^ < 0. In this case, if 
we call the future marginal surface is outer. The future marginal 
surface with Lrfi'^^^ > is called inner. 

The past marginal surface is defined by 6*^"-' — 0, d'^^^ > 0. Similarly, the past 
marginal surface with £^0^"^ > is called outer, and the case with ££6'(") < 
is called inner. 

In some neighborhood of a given codimension-2 surface, we can imaging this 
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codiinension-2 surface belongs to a foliation of this neighborhood such that the 
projection operators are Lie dragging along the null vector This way, we can 
calculate CeO^^^ for an arbitrary normal vector Y. By the similar assumption 
and calculation, we get Ln9^^\ By this consideration, we get the values for 
Ci9^^\ Li,9^"'\ £„0*^^' and LnO^^^ on the given codimension-2 surface (Of course, 
this surface belongs to the different foliations of some neighborhoods according 
to the different deformation vector X's). Actually, similar to get eas. ([65|) and 
([S5)) . it's easy to find the (cross) focusing equations: 

n — 2 
n — 2 

/:„6»(^) = Kr.e^'^^ + uj^Lo" - D^uj'' + ^bi'^n'' -^R- 6l(^)6'(") , 

^ + + 4^ _,(^),("). (94) 

So, for the future marginal surface, the classification of the outer and inner is 
determined by 

Cn9(''> = LJ^LJ' - D,LJ' + ^ad'^n'' - (95) 

on the future marginal surface, while for the past marginal surface, this kind of 
classification is determined by the value 

£<>6l(") = Wcw" + D.w'^ + '^abtn^ - (96) 

on the past marginal surface. Assuming the marginal surface is closed, then, for 
the future outer marginal surface and the past inner marginal surface, above re- 
lations give strong constraints on the scalar curvature R of the marginal surface. 
For example, J eqR should be positive if some energy condition is imposed. 

The so called trapping horizon is the closure of a hypersurface foliated by 
the marginal surfaces. The classification of the trapping horizon inherits from 
the classification of the marginal surfaces fS] [TUJ [TT] . 

So we can imagine that the trapping horizon ^ is a hypersurface which is 
foliated by a family of (n— 2)-dimensional marginal surfaces Sr, where the t & R 
is called foliation parameter of the trapping horizon. Assume X is the so called 
"evolution" vector, i.e., the vector which is tangent to % and normal to Sr and 
satisfies Cxt = 1. The details can be found in, for example, reference j38] . 
Since X is normal to the codimension-2 surface Sr, it can be expressed as 
Xa — Aia — Bua as bcforc. Of course, if we require that the projection operator 
is Lie dragging along X, then X is also a deformation vector. In this case, all 
the formula in previous sections can be used for this X. In following discussion, 
we always assume that the evolution vector X is also a deformation vector. 

Here, it should be noted that the foliation near some marginal surface 5*,- 
has been given in advance. For instance, the neighborhood of Sr characterized 
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by the foliation parameter from t — At to r + Ar is given a priori (and the 
fohation structure of this neighborhood is also provided simultaneously). By 
requiring that qj" is Lie dragging along X, we get the value of LxO^^^ on Sr- 
For another vector, for example, we have no a priori neighborhood of St and 
the associated foliation structure. However, we can get the value of Ce9^^^ on 
Sr by requiring that qj' is Lie dragging along i although the corresponding 
neighborhood and associated foliation structure is not provided explicitly. This 
point has been discussed at the end of Sec. 3. 

For the evolution vector X, some discussions are listed in order: 

(i) . For the reparameter or the relabeling of the foliation t — > t'{t), the 
evolution vector changes a,s X ^ X' = {dr' /dT)~^X. As we have pointed 
out before: According to eg. ([60]) . the relation between CxO^"^^ and Lx'O^^^ 
is complicated due to the existence of the term Y'^DaD°'Xc- So generally we 
have no Cx'0'^^'> = fCxd^^'> if X' = fX (i.e., eq.^). However, since that 
T and t' — t'{t) are both constants on the marginal surfaces, we really have 
Cx'O''^^ = {dr' / dT)^-'- CxO'^^\ Thus, this kind of reparameter does not effect 
our discussion (the classification of the trapping horizons). 

(ii) . It's easy to find that the position of the marginal surface is independent 
of the rescaling of the null frame {£, n} {Xi, n/X} with some positive regular 
function A. However, the value of {CiO'^"^) on the future (past) trapping 
horizon really depends on the rescaling of the null frame. So the classification 
of the outer and inner of the trapping horizons is complicated. However, in the 
case with enough symmetries, the classification of outer of inner of the trapping 
horizon is independent of the rescaling of the null frame. 

(iii) . One can not naively get Cxd^^^ just from the linear combination of 
Ci9^^^ and CnO^^'' listed in eq. (|94|) . The reason is very clear according to eg. ([62|) . 
The values of CxS'^^'^ and Cxd^"''^ are given in eqs. dHSI) and ([55)) . Certainly, for 
some special case (for example, the case of spherically symmetric spacetime), 
CxO'^^^ is really a linear combination of CiO^^'' and £„0*-^-'. 

(iv) . Since the trapping horizon is a hypersurface foliated by the marginal 
surfaces, generally we have 6''^^) = on this hypersurface. Further, since that the 
"evolution" vector X is tangent to this hypersurface, then we expect Cx9^^^ = 
on this hypersurface. From eq. ([S5|) . we get a complicated differential equation of 
the function B. So it's difhcult to get the explicit expression of X. However, in 
some cases with enough symmetries (for example, spherically symmetric space- 
time), the problem becomes very simple. We can find that X can be fixed up to 
a function which does not depend on the points of the codimension-2 surface. 
This will become clear in the FLRW universe. 

6 Horizon Dynamics with Quasilocal Energy 

To study the dynamics of the horizon, we can firstly assume some quasilocal 
energy, i.e., the energy inside the closed codimension-2 surfaces, and then study 
the deformation of this energy. Generally speaking, the deformation of the 
energy can be reduced into a form like the first law of thermodynamics if this 
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quasilocal energy is appropriately selected. Need not to say, this discussion 
heavily depends on the definition of the quasilocal energy. 

6.1 General Spherically Symmetric Cases 

The results in previous sections are valid for any spacetime. In this subsection, 
as an example, we will discuss the spherically symmetric spacetime. Here, in 
this simple case, we will not introduce the quasilocal energy at first. Instead, we 
show how the deformation equation of the expansion can be transformed into 
an interest form where the quasilocal energy will appear naturally. Generally, 
the metric of an n-dimensional spherically symmetric spacetime can be written 
as 

g = P^,{y)dy^'dv^ + r{yf-(,,{z)dz'dz^ , (97) 

where "fijdz^dz^ is the standard metric of an (n — 2)-dimension sphere of radius 
one. Assuming the codimension-2 surface is just the (n — 2)-sphere, then we 
have 

habdx'^dx'' = P^^dyf'dy'' , qabdx'^dx'' = r^j.^dz'dz^ . (98) 
With these identifications, it's easy to find 

1 1 (n — 2) 

Kabc = — gafcVcr = -QabKc , - ^V^r . (99) 

r n — 2 r 

Assuming X and Y are two normal vectors, so, from qabX^ — (labY^ = 0, we get 
that Xadx"' — X^dyi^' and Yadx"^ — Yf^dy^, where Xf^ and are components of 
X and Y in coordinate bases. Further, we require that X satisfies CxQa = 0- 
A short calculation shows that only depends on the coordinates y'^. Thus, 
according to the definition of Da, it's easy to find 

DaXb = q^K'V.Xd = . (100) 

So the two order derivative term of X in ea. (j60l) is vanishing, and the behavior 
of the usual Lie derivative of X in eq. ((62|) is restored. This means that eq. (|60| 
only depends on the direction of X (but not the norm of X). Now, we have 
CxO^^^ = X^Vad^^^ and considering 9^^') = -K^Y", we get 

Cxe'-''^ - X'^V.Y^ ( "^VcA + ^y^X^V.Vcr - Izlx'^Y'VerVcr. 
\ r J r 

(101) 

Substituting above result and ea. (l99l) into ea.(l60|. we have 

(^) X^Y^VaVbr - ^-^^-^^^ [1 - iX^Y,) 

= -^afc [X'^Y^ - h"^ {X^Ye)] . (102) 

Here, we have substituted the value of the scalar curvature of the sphere: 
R— (n — 2) (n — 3) /r^ . Now, let's choose X just be the extrinsic curvature vec- 
tor (or mean curvature vector) of the sphere, i.e., Xa — Ka = —{n — 2)\I a'i'/f- 



24 



Obviously, this selection satisfies the condition that the components of the de- 
formation vector X are only functions of the coordinates y. Considering 

V6(VVVar) = VWbVar + (V6Var)V"r = 2VW6Var , (103) 

and rearranging the terms in eq. (jl02p . we get 

^^^r^Vb [r"-3 (1 - V.rVV)] = %6 [VVF^ - /i''^ (VVFe)] • (104) 
After integrating this equation on the codimension-2 surface, we get 
{n - 2)r2„_2 



(1 - VarVV) 



167rG 

(105) 

where £/ = r2„_2''" ^ is the area of a sphere of radius r, and J?^;, is the energy 
momentum tensor of the spacetime. Here, we have used the Einstein equation 
^ab = SirGXh- By defining 

^ ^ ("- 2)a.-2 ^„_3 (1 _ y^^ya^^ (log) 

IdttG 

and 

V-a = =^afc VV + wVar , W = -^h"'' Xb , (107) 

we can put eq. (|105p into a very familiar form: 

Cy^^ ^ J^i^aY" +wCYy , (108) 

where the normal vector Y is arbitrary, is the volume inside the (n — 2)- 
sphere which is given hy y — n,n-2r"~'^ / {n — 1). It should be noted: we have 
selected a special X which is proportional to the mean curvature vector of the 
(n — 2)-sphere. So we are considering a special foliation of some neighborhood 
of the codimension-2 surface. By this selection of the deformation vector, the 
deformation equation (1601) can be transformed into above simple form. 

The quantity actually is just the so called Misner-Sharp energy insider 
the (n — 2)-sphere ^6j. While V'a is a kind of energy flux which is called energy 
supply. The scalar w is an energy density, and wCy Y in eq. (jlOSp is a work term. 
Ea. (jl08l) is called unified first law, which is firstly found by Hay ward in GR [HI 
[TUl [TTj . Certainly, one can get ea. (|108|) directly from the Einstein equation. 
Here, we have deduced it from the deformation equation of the codimension-2 
surface. 

Above discussion is independent of the null frame. If wc introduce the null 
frame {£,n}, then the mean curvature vector can be expressed as 



if„ = 0Wn„ + 0(")4--^^^Var , (109) 

r 

and other quantities can also be expressed by using the null frame, for example, 
the Misner-Sharp energy S" : 



^ ^ (n- 2)»„_2 ^n-3 

IGttG 



1 + .A^0W0in) 

[n — 2Y 



(110) 
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Certainly, it's also easy to get the expression for the energy supply in this null 
frame. 

Everything is simple when the spherical symmetry exists. However, for the 
general case, the situation becomes complicated. Maybe, the most difficult prob- 
lem is: how to select an appropriate quasilocal energy. In the four dimension, 
to get an equation which is similar to the unified first law (jl08l) , Hayward et. al. 
have used the so called "Hawking mass" [28l [29l |30l [31] . In the next subsection, 
we will give a similar discussion in the higher dimension. However, firstly, we 
have to generalize the Hawking mass (energy) to the higher dimension (n > 4). 

6.2 More General Cases in Higher Dimensions 

Eg. ljllOp depends on the null frame. However, we can put it into another form 



where R is the scalar curvature of the codimension-2 surface. This equation 
still depends on the coordinate function r. By using the area, / e,, of the closed 
codimension-2 surface, we can transform it into a more general form: 



Now, (f does not depend on any applied structure of the codimension-2 surface. 
In static case, a similar mass function has been studied in reference [47] by using 
the (n — 1) + 1 decomposition of the Einstein theory in n-dimension. The energy 
expression (jll2p is interesting: 

(i) . In the four dimension, this energy is nothing but the so called Hawking 
mass (energy) fF7\ . 

(ii) . In the general spherically symmetric cases, it reduces to the Misner- 
Sharp energy in the higher dimension (The generalized Minsner-Sharp energy 
in general Lovelock gravity theory can be found in [FH 15^). 

Tl — 4 

(iii) . The term (/ eqR)/{J e^)"^ has close relation to Yamabe invariant 
(see the chapter 4 of [53] for the details of the Yamabe invariant). For n = 4, 
this term is just Euler number, and its variation vanishes. On the other hand, 
for n > 4, an arbitrary variation of this term also vanishes if the codimension-2 
surface is a compact Einstein manifold without boundary. This can be found as 
follows: For example, we just take the variation to be Cx, i.e., the deformation 
operator with a deformation vector X, then, it's easy to find 




(111) 




-3 
-2 



(112) 




(113) 
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To simplify the expression, let's introduce ^ = J eg. Assuming the codimension- 
2 surface is closed, we get 

C. (h^)=^('-±)^-^fJH('-^. K'X^ \. (U4) 



1-2 



71-2 



where Lxs^ = Cx J — — J ^qiK'^X^). Thus, there are three ways to set 
above variation to be vanished: 

(a) . Obviously, the right hand of above equation is identically vanished in the 

Tl — 4 

four dimension because that J eqR/{J Cq)"^ is just the Euler number of 
some two dimensional closed manifold. 

(b) . If the codimension-2 surface is assumed to be a closed Einstein manifold 

(i? is a constant), then we also get vanished variation. 

(c) . Selecting a special deformation vector X such that K'^Xa is a constant on 

the codimension-2 surface, then we have Lxs^ I + K^X^ = 0. 

To simplify the discussion, here, we introduce a quantity 

^= -l^^^l- (115) 

167rG(f}„_2)~ (n-3) \(/e,)^ 

Once one of above three conditions is satisfied, we have Cx'^ = 0. Of course, 
generally, Cx-^^ 7^ when n > 4. 

From the definition (|112p . the deformation of this energy along the normal 
vector X is given by 

To get this equation, the requirement Cx-^ = is not necessary. However, 
when Hx'^ — 0, we can transform above equation into another useful form 

llxS = { { ^\ Cx^ + ^"^Cx i^-A. (117) 



n-2 



After inserting J^^, the term inside the last bracket in above equation is propor- 
tional to K'^Kc- 

To find V'a and w like quantities as in the spherically symmetric case, we 
need further calculations of the deformation of this energy. Firstly, it's easy to 
find 

-£x {K'K,) = Cxe^""^ 

-K'bcb'Xe - Kc (X^^eK^) , (118) 
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then, we get 



ri — 1 



2 ^ 



R 



KrK' 



[KeX' 



(119) 



where we have used the fact that Cx {K^K^) — 2Kc {X'^VgK'^) and the defini- 
tion of C^f^'^ in eq.©. Substituting this result into ea. pi6l) , after some algebraic 
calculations, we get an important formula of this section: 



1 




K,X^ 
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K^DcD^Xe 



i {'^abh'^') ■ (K^X^) i , 



(120) 



where L = £/ 



£ ^ 



and 



L 



16nG{n - 3) 



R- 



n — 3 
n - 2 



KcK' 



(121) 



Here, f is a quantity like an energy density. In fact, from the definition of the 
energy S", we have = / CqS- Ea. (ll20p gives the deformation of the energy ^ 
inside the closed codimension -2 surface. Some remarks are listed in order: 



(i). 



(ii). 



If n = 4, the energy (a is the Hawking mass, and ea. (|120l) reduces to the 
one given by Bray et.al [31]. Of course, we can also consider the cases with 
cosmological constant as in |31| . According to previous discussion, in this 
case, we have CxJ^ = because .J€ is actually a topological quantity 
now. However, we have to point out: The four dimension is not so special 
if we put the deformation of the energy into the form (I120p . 

In the general spherically symmetric cases, C^^ and K'^DcD^X^ are van- 
ishing. Further, the first term in the right hand of ea. (|120p is also van- 
ished, then ea. (|120p reduces to the one given in ea. dlOSp . This is an 
expectable result. 



(iii). It should be noted that the vector X can be changed arbitrarily in this 
equation. So, if we select X such that K'^X^ = 0, then we get 



Cx<^ 




K^DrD^Xe 
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(122) 



Further, when X is just the dual vector of Kc, then, K'^Xc automati- 
cally vanishes and the term K^DcD'^Xe is a total derivative on the closed 
codimension-2 surface (see ea.([7T|)). So, in this case, the term K'^DcD'^Xe 
can be omitted. 

(iv). It's easy to find: the first term in the right hand of eq. (|120p vanishes if 
f is a constant on the codimension-2 surface. This means, on the equi-f 
surface, we always have 



However, unfortunately, the marginal surface (on which we have KcK'^ = 
0) is generally not an cqui-£ surface unless it is an Einstein manifold. Of 
course, we also have Hx'^ = if the codimension-2 surface is an Einstein 
manifold. 

(v) . We can get the result ()123p by another method: Tuning the vector X 

such that KcX'^ is always a constant on the codimension-2 surface. By 
this selection, we have CxJ(^ — 0. The requirement in the item (iii) 
{KcX'^ = 0) is just a special case of this method. 

(vi) . It's interest to study the monotonicity of this energy ^ as in [3TJ |47] . 

Although J(f is not a topological term for n > 4, the deformation of the 
energy (a, i.e., eq. (|120l) . really has a very similar structure as in the four 
dimension. So, it's possible to get some monotonicity behavior of this 
energy. However, this is far beyond the aims of this paper, and we will 
not give further discussions here. 

Above discussions tell us: The situation becomes a little bit complicated 
when the codimension-2 surfaces are not Einstein manifolds. In the following 
discussion, to simplify the problem, we will firstly discuss the case in which the 
Einstein condition is imposed. Since the (n — 2)-dimension Einstein manifolds 
are always constant curvature spaces when n < 6, so, in lower dimensions, the 
generalized Hawking mass (|112|) has no big difference from the Misner-Sharp 
Energy in the spherically symmetric cases. However, the problem becomes 
interesting when n > 6 because in these cases the Einstein manifolds may be 
inhomogeneous. To discuss the situation without the Einstein condition on the 
codimension-2 surface, we have to consider some special deformation vector X. 
We will study this situation at the end of the next subsection. 






(123) 
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6.3 Dynamics of Trapping Horizon 



The discussion in previous two subsections have given enough preliminaries to 
study the dynamics of the trapping horizon. Firstly, let's consider the general 
spherically symmetric cases. 



6.3.1 Spacetime with Spherical Symmetry 

In previous sections, with the spherical symmetry, we have shown that the 
deformation equation can be transformed into the form of eq. (llOSp . By choosing 
X (This corresponds to the selection of Y in eq. (|108p ) to be the tangent vector 
of the trapping horizon on which VarVr = 0. More precisely, we choose X to 
be the evolution vector of the trapping horizon. Then, on the trapping horizon, 
we have £x(VarV°r) = 0, so ea. (jl05l) becomes 

Cx'^=[^)CxS, (124) 

where 

K (n — 3) c 

2^^^^' ^"^AG- ^^^^^ 
This K is just the "effective surface gravity" studied in [6l[61]. Another inter- 
esting surface gravity is defined by 

^V-a^'^ (126) 

which has a form 



K _ AG 
2tt ^ n~1 



(127) 



In above expression of the surface gravity k, the energy S takes value on the 
trapping horizon. In the four dimension, this surface gravity is just the one 
defined by Hay ward |^. By this k, the evolution of S on the trapping horizon 
becomes 

CxS=[-^)CxS + wCxr ■ (128) 

This is a first law like equation. Obviously, k is a constant on the marginal 
surface {DaH = 0). However, generally, it's not a constant on the full trapping 
horizon. The reason is: generally, we have no Cxi^ = 0. So k may evolve on the 
trapping horizon. 



6.3.2 General Cases with Einstein Condition 

In the more general cases where the codimension-2 surface is Einstein (so Hx-^ = 
0), we have eqs. ()117p . (|120|) and (|123|) which have been found in previous sec- 
tions. Ea. (|117p is simple. If we set X to be the evolution vector of the trapping 
horizon, then we get 

LxS^^LxS, (129) 
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where the "effective surface gravity" k and the entropy S are defined as 

'?i-3\ / S" 



27r 



n 



* ~ 4G 



(130) 



Of course, this kind of surface gravity is effective, and usuahy it can not reduce to 
the surface gravity of the corresponding stationary spacetime. By this definition, 
it's easy to find [n - ={n- 2)fS with f = ii/2TT |T]. 

For the marginal surface, we have K'^Kc — 0, so ea. (ll20p reduces to eq. (|123p . 
This is a kind of energy balance equation. Actually, from the Einstein equation 
%b = 8nG3^ab, eq. (fT23| becomes 



"2 i^abh'"'' 




(131) 



In the general spherically symmetric cases, in above equation, the terms re- 
lated to the energy-momentum tensor 3^ab give A^aX"- -\- wCx'f ■ However, 
in the general cases, there are two additional terms corresponding to C^^ and 
K^DcD'^Xe- Of course, without enough symmetry, it's also not easy to complete 
the integral, and we can only write the contribution of the energymomentum 
tensor .!^ab into an integral form 



(132) 



where w'^"^^ and i/j^™) are defined as 



w^"'^ ^ --Xbh''\ ^^abH' + W^-'^^Ha. (133) 



Obviously, w^™) is the same as the one in the spherically symmetric case. The 
vector H"" plays the role of V"r in the spherically symmetric case, which is 
defined ail 

Ha = , (134) 



Eq. (|132p conies from the matter fields. While the terms represented by C^f^'^ and 
K'^DcD'^Xe can be understood as the contribution from the change of a gravi- 
tational field. In some sense, we can understand them to be some gravitational 
radiation. It's obvious that shear tensor C^j'' provides an energy supply 



SttG 



C/^cd 
cda'^ b 



2 ^ 



(135) 



■^By using the normal vector H", we can also define a generalized Kodama vector as: 
^abH^ 1481 . In the spherically symmetric case, it's possible to define the surface gravity from 
this vector, for example, see reference 1521 . 
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as the usual energy- momentum tensor ^ab- However, it has no contribution to 
w. Instead, it seems that the term K'^DcD'^Xe plays the role of w. Actually, if 
K'^Xf^ is nonvanishing, formally we can define 



,(9) - 



1 



SttG 



(136) 



With these identifications, the contribution of the gravitational radiation is 

j(^^(9) +uj^a)Ha)x\ (137) 
So the deformation of the energy S" becomes 

Cx'^ = / { (V^i"^ + ^i'O + ("'^"^ + ""^'0 ^4 • ^^^^^ 

Similar to eq. (|126p . we hope the energy supply could provide a similar definition 
of some surface gravity as the one in eq. ()127p . Assuming X to be the evolution 
vector of the trapping horizon, and using eq. ()129p . we get 



= - / e 



(i^aX") 



(139) 



Unlike the case with the spherical symmetry, usually, we can not get a surface 
gravity which is a constant on the marginal surface. Actually, by defining 



K 

2^ 



4G 



,„-3)(^)-.( 



,,(9) 



and considering Cx^q = d'^-^^eq — —{K°'Xa)€q, at most we have 
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^xs , 



(140) 



(141) 



where s = Cq/AG can be simply understood as the entropy associated with 
the area element of the marginal surface. Obviously, the surface gravity (|140p 
reduces to the one in eq. (|127l) if the spherical symmetry is restored. So ea. (|138l) 
becomes 



Cx'^ 



wCxv . 



(142) 



where w = w'^"^^ +w^^^ and Cxv = O^^^eqL/ (n— 2). In the spherically symmetric 
case, above equation reduces to eq. (|128p . 

Above discussions are independent of any local frame. It's also clear that 
above description does not depend on the relabeling of the foliation of the trap- 
ping horizon {X — > fX with some relabeling factor /). Since only KcK'^ = 
and Cx{KcK'^) — (on the trapping horizon) have been used, so all discussions 
are valid on any kind of trapping horizon (future or past, outer or inner). 
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6.3.3 Expression in Null Frame 

However, to make the problem easy to understand, it's better to introduce some 
null frame. For an arbitrary null frame {i,n} with f = — 1, we can express 
X and K to be X° = At' - Bn'' and K"" = 6'(^)n'^ + 6i(")^° respectively. It's 
easy to find K^'Xa = -A6'(^) + B9^"\ From cq.§TI, we get 

+ e*'") {D^D^B ~ 2uj''D,,B - BD^uj'^ + Buj^uj") . (143) 

So, on the future trapping horizon [B^^'i — and 6*^"'' < 0), when K'^Xa ^ 0, we 
have 



K'^D D'^X 1 

^ — (D'D.B - 2uj''D^B - BD^Lo" + Buj^uj") 

K'^Xa B^ ^ c ^ c ) 

Then, ea. (|136l) becomes 



(144) 



w 



(a) 



1 



4 ) {D^DcB - 2uj''DcB - BDcOj" + BuJcUj'') . (145) 



8ttG \B^ 

In the null frame, the energy supply (I135p is also very simple: 



1 



inG \n-2 

Similarly, the energy supply of the matter fields becomes 

L 



n-2 



(146) 



(147) 



After substituting above expressions into ea. (|138p . we get 



+B^abtn'' 



L 



n-2 
1 



A [ ^abt& + 
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{D^DcB - 2lu'DcB - BDcUj'^ + BlJcLj'' 



.(148) 



This result is valid even in the case where K'^Xa — 0. If K'^Xa = 0, we can't 
define w'^^^ as eq. ()136p . Actually, now, the deformation of the energy reduces to 
eq. (|122p . Further, since K°-Xa = BO^"^^ on the future trapping horizon, so, to 
require that K°'Xa = 0, we have to set X° = At (i.e., the trapping horizon is 
null). Thus, eg. p^ implies that K'^DcD'^Xe = 0. Therefore, in the nuU frame 
which just has been used, eq. (|122[) becomes 



L 



n-2 



A Xbtt + 



(149) 
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Obviously, we can get this result just by setting i? = in eq. (|148p . 

Generally, the evolution vector X has to satisfy CxO^^'' = on the future 
trapping horizon. Then, considering ea. (l55|) and K'^Xa — B9''^'^\ the relation 

(fT48l) becomes Cx<S = (^) {^) Cxs^- This is nothing but eq. (|117l) taking 
value on the trapping horizon. This resuh imphes that the energy S does not 
evolve on the null trapping horizon on which we have K'^Xa = S6i(") 0. This 
point can also be directly found from eas. (|55|) and (jl49l) . 

Above discussion is valid for any null frame {£, n\ which satisfies the relation 
£°na = —1- Sometime, one can rescale the null frame {i,n\ such that 0^"^ 
satisfies HOI 121 121] 

= -^^-^ < 0. (150) 
ij 

It's easy to find this requirement also means = — 1 on the marginal surface 
of the trapping horizon. Now, eq. (|148p is very simple: 



1 



SttG 



+ B .%b(.°-rf + 



SttG 



(151) 



where Cc = Wc — Dc\^B. Since the codimension-2 surface is closed, all possible 
total derivatives can be omitted. So after setting ^^^"^ — — (n — 2)/L, we can 
omit the total derivative terms in eq. ()148p . By this consideration, eq. (|149p also 
becomes simple: 



CxS^ 



SttG 



(152) 



However, it should be noted here: After selecting this kind of null frame, we 
can not get result (|152p just by setting B = in ea. (|151l) . This is because there 
are two Da \nB terms in cq. (|15ip . 

In this null frame, cq. (|15ip tells us that the evolution of the energy S on the 
future trapping horizon is decomposed into two pars: 

(i) . The contribution of the usual matter fields — ^6^°^" and ^b^°rt° ; 

(ii) . The contribution of the gravitational radiation — '^^^^'^^ ^^"^ CaC- 

We will give more discussions on the gravitational radiation at the end of this 
section. For the past trapping horizon, it's also easy to get similar results like 
eas. (|15ip and (|152p . Since the procedure is similar, we will not give further 
discussions here. 



6.3.4 Codimension-2 Surface without Einstein Condition 

Without the Einstein condition, the problem becomes complicated even in the 
four dimension. For an arbitrary null frame {i, n} with tan"" = — 1, the evolution 
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of the energy becomes 
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{D^D^B - 2uj^DaB - BD^Lo" + Boj^oj") 



(153) 



Generally, the evolution vector X°- = Ai°'~Bn°' does not satisfy the requirement 
that K^Xc = -B6'(") is a constant on the marginal surface. However, in the case 
of i? > 0, we can always rescale the null frame such that 



B6i("' = -- 



L 



(154) 



is a constant on the marginal surface (Of course, L is a constant on the marginal 
surface). It's easy to find this requirement also means CxL = 1 (This selection 
is also used in [29l[30]). Thus, by using this special null frame, the evolution 
of the energy becomes 



Cx<^ = 




8^G '^^'^ 



b , CrX" 



8ttG 



(155) 



where C,c is the same as the one in ea. (ll51|) . Considering that X has to satisfy 
eq. (p5)) with 6l(^) = and CxO^'^'^ = 0, we have 
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CxL 



eqR = 



n — 3 
n-2 



S 



(156) 



where S takes value on the marginal surface. To get above result, we have 
inserted the relation CxL — 1. It should be noted here: when K'^Xa is a 
constant on the codimension-2 surface, then we have CxJ^ = 0, so eq. ()116p 
can be transformed into eq. (|117p . Then, on the trapping horizon, the second 
term in the right hand of ca. (|117p vanishes, and then we get above equation. 
Eq. (|155p also shows the evolution of the energy inside the marginal surface can 
also be decomposed into two parts: The contribution of the matter and the 
contribution of the gravitational radiation . 

When the tapping horizon is nuU, we have K'^X^ = BO^'"-^ = 0. Eq. p351) 
reduces to ea. (|149l) . So, on the null trapping horizon, the energy inside the 
marginal surface does not evolve regardless the marginal surface is Einstein or 
not. 

Here, some remarks are listed in order: 
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(i) . The deformation of the generahzed Hawking mass (I112p . i.e., eq. (|120p . 
has a very similar form as the one in the four dimension. When the marginal 
surface is an Einstein manifold, the evolution of the generalized Hawking mass 
(|112p on the trapping horizon is given by eq. (|129p . This result can be decom- 
posed into the form (|142p . In the null frame, this result can also be expressed 
in eas. (|15ip and (|152p . In the general case without the Einstein condition of 
the marginal surface, by choosing some special null frame, the evolution of the 
generalized Hawking mass can also be decomposed into the matter field part 
and gravitational radiation part. 

(ii) . The gravitational radiation is carried by cr^^'' (or a^^^ in the past case). 
The number of the degrees of freedom of this tensor is n{n — 3)/2. This is just 
the number of graviton polarizations in n-dimension. In the previous discussion, 
we think that the term which corresponds to CaC" is also a kind of gravitational 
radiation. However, we have to point out: the detailed physical meaning of the 
term CaC" is still unclear. Certainly, this term really comes from the change of 
the gravitational field, so it's reasonable to regard it as a term of gravitational 
radiation. 

(iii) . The problem of the angular momentum: The deformation of the Hawk- 
ing mass or it's generalized version in ea. (lll2p in some sense are not sufhcient 
to describe the full dynamics of the trapping horizon. One has to consider the 
angular momentum separately j541 155| . We can define the angular momentum 
to he = J eq{(f)°-uja) for some tangent vector (j) which satisfies Cx^'^ — and 
Dqc/)" = 0, then, the deformation of the angular momentum can be studied by 
eg. ([5^ (or ([M]) '). In the four dimension, there is only one possible angular mo- 
mentum. However, in the higher dimension, the situation is complicated: It's 
possible that there are several (not single) angular momentums associated with 
the trapping horizon. 

(iv) . Generally the surface gravity defined in eq. (|140p is not a constant on 
the marginal surface, i.e., DaK, ^ 0. In the spherically symmetric case, this 
kind of surface gravity is really a constant on the marginal surface. However, 
generally, it still evolves on the trapping horizon, i.e., Cxn 7^ 0. This is very 
different from the static case in which the surface gravity is a constant on the 
horizon. Although we can get some energy balance like equation (for example, 
eq. (|15ip ). however, frankly speaking, the definition of the surface gravity of 
general nonnuU trapping horizon is still an open problem |49| . Of course, we 
hope the definition of the surface gravity can give some physically acceptable 
result when the system is almost equilibrium. 

(v) . For the null trapping horizon, the generalized Hawking mass does not 
evolve on the trapping horizon. Further, if null energy condition is imposed, 
all terms in the right hand of eq. (|152p have to be vanished. So there are no 
dynamical version first laws associated with the null trapping horizons. To 
study the dynamics of the null trapping horizon, we have to consider other 
method, for example, the phase space version first law |13[ [Ml [TSl ITBl [17] . 
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7 Horizon Dynamics without Quasilocal Energy 



Without selecting some quasilocal energy inside the codimension-2 surface, we 
can also study the dynamics of the trapping horizon. This kind of discussion 
heavily depends on the deformation equations of the expansions and the SO{l, 1) 
connection we have get in section 4. Eas. ([75)) and ([M]) are key equations to 
construct some first law like equation in this formalism. Certainly, one of the 
most important problems in this formalism is also the definition of the surface 
gravity. As mentioned at the end of last section, to make the problem easy 
to understand, it's better to discuss some near equilibrium state at first. So, 
in this section, we discuss the slowly evolving trapping horizon proposed by 
Booth et.al. |20l EU [22l |46] . This theory just describes the trapping horizon 
which corresponds to the near equilibrium state. We will generalize the slowly 
evolving future outer trapping horizon to the past trapping horizons (sometime 
there are no future trapping horizons in a given spacetime). This is important 
to study the FLRW universe because the past trapping horizon is important in 
this case. 

To make the discussion clear, here, we put the deformation equations of 
the expansions (After selecting X to be the evolution vector, these are actually 
evolution equations), i.e., eas. (j65p and (1681) into simple forms: 

C^e^'^^ = -D^D^B + 2u''DrB - BC^e'^^^ + ACiO'^'^^ , 

/:^6'(") = DcD-'A + 2u;'D^A + ACe9^"^ ~ BLn9^''^ . (157) 

For the future trapping horizon, we have CxO^^^ = 0, so the relation between 
A and B is given by a two order differential equation of B. Similarly, for the 
past trapping horizon, we have CxO^"''^ = 0, and the relation between A and B 
is encoded in a two order differential equation of A. 

7.1 Equilibrium State 

The theory of the so called slowly evolving horizon is proposed to describe 
the dynamics of some horizon which corresponds to the near equilibrium state. 
So, to study the slowly evolving horizon, firstly we have to study the horizon 
which corresponds to an equilibrium state. Actually, in black hole theory, the 
equilibrium state is described by the dynamical behavior of some horizon which 
is null. Of course, this means the vector X is null. For the null future trapping 
horizon, we can assume Xa = Aia, and then the evolution equations of the 
expansions become 

CxO'^'^^ = ACee^^'^ , CxO'^"'^ = D^D'^A + 2uj''DcA + ALfi^^"^ . (158) 

Similarly, for the null past trapping horizon, by setting Xa = —Bua, the evolu- 
tion equations are simplified to be 

/:^6'(") = -B/:„6'(") , £xO^'^'> = -D^D'B + 2uj'D^B - BCnO^^^^ . (159) 
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Firstly, assuming the null energy condition is satisfied, then by using Cxd^^^ = 
and the (cross) focusing equations (I94p . on the null future trapping horizon, 
we have 

frit^^O, ^ad-'i'^O, (160) 
and on the null past trapping horizon, we have 

= ' ^afen-^n" = . (161) 

Above equations tell us: K^^^ = on the null future trapping horizon, and 
^ab = on the null past trapping horizon. Further, 'Sah^"'^ = and %bn°rt^ = 
just imply that there are no matter flux across the codimension-2 surface. 

Secondly, if we also require that q^^ch^ = on the null future trapping 
horizon and q^'Scb'f^ = on the null past trapping horizon, then, from the Co- 
dazzi equations ([M)) and relation (1551) . we get qa q^'^ (.'^'toebcd = on the null 
future tapping horizon and q^ q^'^n'^^ehcd = on the null past trapping horizon 
{q^^cb^^ = or q^'i^cb^^ = can be satisfied if dominant energy condition is 
assumed.). 

Finally, from above requirements and cq. (|80p . it's easy to find 

LxUJa - DaKx = (162) 

on the null trapping horizon. This result does not depend on the selection of A 
(B) in Xa = A£a ( Xa = —Bua)- According to ea.ipn]). it's also independent 
of the rescaling of the null frame {£,n}. Additionally, if one requires that uJa 
does not evolve, i.e., Cx^a = 0, then, from above equation or ea.(|Ml). one 
gets DaKx = on the codimension-2 surface for both cases (future and past). 
This means kx is a constant on the codimension-2 surface. Furthermore, if 
Cx^x = is required, then kx is a constant on the null trapping horizon 
(regardless future or past). 

In fact, the requirement that kx is a constant on the null trapping horizon 
gives some constraints on the function A (B). In the case of future, this can be 
found from following equations: 

DaUX = KiDaA + ADaKe = , 

CxKx = A{Ke£iA + ACeKi) =0. (163) 
Similarly, for the past null trapping horizon, B has to satisfy 

DaKx = -KnDaB - BDaKn = , 

CXKX = B{K,,£nB + BCnKr,)=0. (164) 

Since, until now, the null frame {£, n} can be arbitrarily rescaled: {£, n} — > 
{Xi, n/A}, generally, ni (k„) is not a constant on the future (past) null trapping 
horizon. Given a null frame {£,n}, we can always find some A (B) to satisfy 
eq. (|163l) (([164])). Obviously, this A (B) is not unique. Inversely, given an 
evolution vector X", we can always find some null frame {£, n} (also not unique) 
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to satisfy eas. (|163p and (I164p . So eas. (ll63p and (|164p give the relations between 
the fohation structure (X) and the nuU frame. 

For example, if we have select some A such that ki (k„) fulfills the require- 
ment that K{ (k„) is a constant on the future (past) null trapping horizon, 
then, from eq. p^ f (ITM)) ). we have to set A (B) such that Da A = LgA = 
{PaB = CnB ~ 0) on the null future (past) trapping horizon (we have assume 
that K(, and both nonvanished). Obviously, these conditions can be sat- 

isfied if we choose that A (B) is a constant. In this case, it's easy to find kx 
satisfies (From eqs.fES]) and (gH), it's easy to find q^X'^V^Xc = 0.) 

X'^V aX^ ^ ±KxX^ , (165) 

where corresponds to the case of future, while " — " corresponds to the case 
of past. This is just the usual formula to define a surface gravity. Obviously, 
this kind of surface gravity is defined up to a constant coefficient. 

Actually, for a given function A (B), we can set A to be proportional to A 
with a constant coefficient (or A which is proportional to l/B), and rescale the 
null frame such that X°'WaX'' = ±kxX'' (with constant kx) is satisfied under 
the resulting null frame. 

It should be noted here, we have not discuss the relabeling of the foliation 
until now. The vector X is rescaled by a factor /(r) = {dr' /dr)^^ if we relabel 
the foliation: r — )■ t'{t). For the null trapping horizons, this just means that A 
or B is rescaled by the factor /. So we can use the same procedure as before to 
find a special null frame such that the relation (|165l) is always held. 

From above discussion, we find: to make that kx is a real surface gravity 
given by eg. U65]} . it's necessary to select a special null frame to match the 
given foliation structure of the null trapping horizon (i.e. X). Certainly, the 
physics the horizon should not depend on the relabeling of the foliation and 
rescaling of the null frame. So different X's (with the preferred {£, n} and the 
corresponding kx) in some sense are physically equivalent. In summary, the 
null trapping horizon can be characterized by an equivalent class which can be 
expressed by a triplet 

X,{l,n},K.x . 

Of course, to foliate the null trapping horizon, the simplest way is select X° = 
(or X"" = — rt°). With this selection, ea. (|165l) is automatically satisfied. 

Above statement has close relation to the isolated horizon (especially in the 
case with constant A [B)) defined by Ashtekar el.al. [13]. In those cases, one 
mainly focuses on the future outer trapping horizon. Here, we also study the 
past trapping horizon. Our discussion is independent of the selecting of the 
metric of the null hypersuface (degenerate). In fact, this is a rough way to 
reconstruct of the isolated horizon [22] • 

Conclusively, on these null trapping horizons, there are no gravitational ra- 
diation and matter fiux, and kjc's are constants. These properties correspond 
to the equilibrium state of the thermodynamics of the horizon. Further, now 
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eas.([75]) and (15^ just mean 

(^^^CxS^O, CxJ^ = 0, (166) 

where S* ~ / and ^ J eq {4>°'U]a) can be explained as the entropy and the 
angular momentum associated with the null trapping horizons. Since there are 
no matter flux and gravitational radiation across the null trapping horizon, these 
physical quantities do not change along the trapping horizons as expected. So 
there are no dynamical version of the first law associated with the null trapping 
horizons. Actually, this point have been found in the Sec. 6.: the energy S does 
not evolve on the null trapping horizon. In fact, one can study the first law of 
the null trapping horizons by using the phase space method ^13l [TH [HI [HI [17] . 

7.2 Near Equilibrium State 

The near equilibrium means that X is almost a null vector. This suggests that 
Xa = Ala — Bua slightly deviates from a null vector. However, to characterize 
this small deviation, it's not sufficient to set one of A or i? to be very small. The 
reason is that there are two ambiguities for the vector X we have mentioned 
several times: 

(i) . There are freedoms to choose the null frame {£, n} -> {A£, n/X] for some 
positive function A, and this makes an ambiguity for B/A or A/B by factors }? 
or 1/A2. 

(ii) . The relabeling of the foliation of the trapping horizon r — > t'{t). 
This makes an ambiguity that X can be rescaled as X X' = fX with 
f = {dT'/dr)-\ 

We can assume that the norm of X (or XaX"^ — 2AB) approaches zero such 
that X is almost a null vector. This will eliminate the ambiguity of rescaling of 
the null frame. However, since X' and X are essentially equivalent to describe 
the trapping horizon, so the norm of the evolution vector will get a factor |/| 
if we use X'. We need some procedure to carefully treat the rescaling of the 
foliation. 

For the null trapping horizons, we can always find some special null frame 
to satisfy the requirement of X^VbX" = zLkxX"" (with a constant kx), and 
above two ambiguities in some sense are fixed to get the preferred null frame 
and the corresponding surface gravity. However, for the general case of non-null 
horizon, the situation is very different. From eq. (j42[) . it's not hard to finc(f| 

X'VbXa = tix {eabX^) - Da (AB) + (CxA) la - {CxB) Ua . (167) 

Obviously, now, it's impossible to get X''Vf,X" = ±kx^" by rescaling the null 
frame. So, principally, the strategy of the rescaling of the null frame in the 
null horizon case is meaningless for the non-null trapping horizon. Actually, 
in this case, we do not know how to select a preferred null frame to define a 
surface gravity associated with X. However, mimicking the null case, for a given 

•^This equation reduces to the result given in p?8j when A = 1. 
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foliation of the trapping horizon, by rescaling the null frame, we can always set 
A (or B) to be a constant for the future case (or the past case). Without losing 
generality, we can set A ~ 1 (or B ~ 1). This procedure is similar to the null 
cases. Thus, for the future trapping horizon, we have 

X" = r - Cn" , (168) 

while for the past trapping horizon, we have 

X" = Cr - . (169) 

It should be noted here, the functions C"s depend on the relabeling of the folia- 
tion of the trapping horizons. Actually, the most general form of the evolution 
vector can be expressed as 

X" = /(r)X° = /(r) (At - Bn") . (170) 

with some null frame {£, n}. One can regard that the factor / comes from 
the relabeling of the foliatior0. For the future case, redefining the null frame 
r -)> ^'^ = {fA)t and = n^/ifA), we can put X into the form ([TM|) 

with 

C = fAB. (171) 

Similar situation happens in the past case. By this, the relabeling of the foliation 
is encoded in the function C. So it seems that one can relabel the foliation of the 
trapping horizons such that C's to be arbitrary nonvanished values (Obviously, 
for null cases, C's are always zero and independent of any foliation structure). 

Since for any kind of trapping horizon, we can always take |C|'s to be arbi- 
trarily small value by selecting the foliation parameter r, so the assumption of 
small |C|'s is still not enough to discuss the slowly expanding behavior of the 
trapping horizons. To describe the almost null property, we need some quantity 
which is independent of the relabeling of the foliation and the rescaling of the 
null frame. We will discuss the cases of future and past separately. 

• For the future trapping horizon. Booth et.el. |20[ [3TJ [HJ give three 
slowly expanding conditions (here we gives a generalized n-dimension version): 

(F-i). The so called evolving parameter e ^ 1 with 



— = max 



\C\ ( ||cr(")f + (87TG)Xbn''n'' + ^— 6»(")6'('' 

\ 71 — 2 



(172) 



*Here, to make the discussion clear, we also require that A (or B) can not be further viewed 
as a relabeling factor. For example, A (or B) is not a constant on the codimension-2 surface, 
so it can not be absorbed into /. Certainly, in some special cases, A and B are both constants 
on the codimension-2 surface. In this case, we can absorb A (or B) into the factor /, and the 
evolution vector is simply assumed to be X" = f{i'^ — C'n°') (or X" = f(C'i'^ — n")) with 
some function C". 
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(F-ii). The Ricci scalar, the SO{l, 1) normal connection and the energy-momentum 
tensor satisfy 



WiUaf and (87rG)^abrn^ ^ ^2 



(F-iii). The derivatives of horizon fields are at most the same order in e as the 
(maximum of the) original fields. For example, 

\\DaC\\^^, \\DaD,C\\^^. 

Here, || • || is the norm of (tangent) tensor fields on the codimension-2 Riemannian 
manifold, while | • | is the absolute value of some scalar. The quantity L is some 
length scale of the codimension-2 surface. For example, the radius of the closed 
(n— 2) manifold: L — {£/ /fi„-2)"^ which has been defined just bellow eg. (|f20p . 
Cm is the maximum value of |C| on the codimension-2 surface. The relation 
E < F means E < k^F for some constant of order one. 

Before using these conditions, we give some discussions: 

Firstly, e defined in the condition (F-i) is independent of the relabeling of 
the foliation and the rescaling of the null frame. This can be easy found from 
the expressions of e in the null frame {^, bellow ea. (|170p : 



— = max 



This is just what we hope to find: The parameter should not depend on the 
relabeling of the foliation and the local frame. By this consideration and the 
definition of e, the requirement of e 1 in condition (F-i) essentially gives some 
constraint on the dynamical behavior of the codimension-2 surface. 

Secondly, the evolution vector X is not arbitrary, and it has to satisfy 
Cx9^^^ = 0. This gives a differential equation of C: 

CxO'^'^^ = -DcD^C + ii^j^D^C - CUB^^^ + dB^^^ = . (173) 

This is just a special case of the first equation in (|157p . Furthermore, considering 
Daj = and the expression of LgQ^^^ in eg. dMl) . the relation LxB^^^ = does 
not give any constraint on the relabeling factor /. This means eq. (|173p is 
actually a equation of AB in (|17ip . Thus, the function C//^ is determined 
by the geometrical behavior (both intrinsic and extrinsic) of the codimension-2 
surface. Remembering that e ^ 1 has some requirement on the geometry of the 
codimension-2 surface in the trapping horizon, so C//^ (or AB) is also required 
to satisfy some condition by the behavior of e. Actually, in some simple case, 
we can find the explicit relation between AB and e. This situation happens in 
the case of the FLRW universe, and we will find it in the next section. 

So the relabeling of the foliation is arbitrary until now. However, if we 
require 

\C\ < <? , (174) 
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then, the function / has to satisfy some condition. This corresponds to some 
special selection of the foliation parameter r, and then we can not relabeling 
the foliation arbitrarily. In the following discussion, we will always assume this 
condition on |C|. 

Now, let's discuss the implication of the conditions (F-i), (F-ii) , (F-iii) and 
(|174p . From these slowly expanding conditions and ea.(|M]). it's not hard to find 
on the future trapping horizon, we have 



(175) 



then, considering eq. (|173|) and |C| ^ e^, one gets that 



Cf, 



STrG^abrr ^ — (176) 



is satisfied on the future trapping horizon. If the null energy condition is as- 
sumed, the two terms in the middle of above equation are both nonnegative. 
It's also easy to find 

Obviously, the first term is the order of e, while the terms proportional to C 
is the order of e^. Since e ^ 1, we have K^^f ^ « a'fjj. Similarly we have 

(Y) (£) 

^ab ~ "ab ' With the same discussion, we get 

XbX'^Y' = Xbi^'i'' - C^Xbu'^n'' = Xbi^'i' + ^(e^) ■ (178) 

Here, the vector Y is the dual of X, i.e., ^ + Cn"^. For other quantities, 
one can also get reasonable approximations. 

In addition to the null energy condition, let's assume that energy -momentum 
tensor also satisfies dominant energy condition. This means, for every future- 
pointing causal vector field Z"', the vector field —^abZ'' must be a future point- 
ing causal vector. This assumption of energy-momentum implies 

g'^'XbZ'^c.dZ'' = Ik.'^fcc^lP - 2{Xbrz'>){£r^dn^Z'') < . (179) 

By selecting = we get 

and then 

ka'^bdl^^- (180) 

From the Codazzi equation ([55]) and (PSj) and above result, we get 

ha^q'^e^'^ebcdW ^ ■ (181) 

For the equilibrium state, the function C is identically vanished. So e is also 
vanished. Obviously, above conditions of the tensors (on the future trapping 
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horizon) give a^^.^ = 0, = 0, q^^ .%cl^ = and q,^ q'^H'^'^ebcd = 0. These 

are just the requirements for the case of the null trapping horizon. 

During above discussions, all the conditions and results are focused on some 
given marginal surface. However, to study the evolution of this marginal surface, 
these conditions are still not enough. Remembering in the case of future null 
trapping horizon, to ensure that some physical quantities (the area and the 
angular momentum associated with the horizon) do not evolve, we have required 
the condition Lx^a = and LxKx = 0. These just mean that uja and kx do 
not evolve respect to the evolution vector Similarly, here there are also 

slowly evolving conditions : 

(F-i'). WCxuJaW and \Cxi^x\ < e/i'; 

(F-ii'). \Cxe^'"^\ ^ e/L2. 

From eq. ((80|) and the restriction condition of the fields (and their derivatives) 
on codimension-2 surface, it's easy to find 

WDanxWdijj^. (182) 

Considering that the absolute value oi Cx^x also satisfies this conditions, then 
Kx is almost a constant on the trapping horizon (at least for some finite interval 
of the foliation parameter r). Thus, kx can be expanded as 

= Ko + , (183) 

where Ko is the leading term of the expansion. 

The condition (F-ii') just requires that 6''^"^ also evolves slowly. It should be 
noted here: generally, in the definition of the future trapping horizon, there are 
no requirements on CxO'''^^ although that the sign of 0*^"^ has close relation to 
the classification of the future trapping horizon. 

• For the past trapping horizon, we can gives similar conditions to describe 
the slowly expanding properties: 

(P-i). The evolving parameter e ^ 1 with 



— = max 



\C\ 



(P-ii). The Ricci scalar, the 50(1, 1) normal connection and the energy-momentum 
tensor satisfy 

\\uj,f and (8^G)^a6rn''^ -1; 

(P-iii) . The derivatives of horizon fields are at most the same order in e as the 
(maximum of the) original fields. For example. 
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Similar to the future case, we can choose some foliation parameter such that 
|C| is small and satisfy the condition |C| ^ e^. 

Substituting these conditions into the expression of ££6''"-' in ea. (|M| . it's 
easy to find on the past trapping horizon, we have 



^ ^ ■ (184) 



then, by setting B = 1 and A — C in eas. (|157p and considering Cxd''"'^ = on 
the past trapping horizon, we get 



S^G^^fen'^n" ^ — . (185) 



Certainly, by assuming the null energy condition, the two terms in the middle 
of above equation are nonnegative separately. It's also easy to find 



and 



K^I^=CKi'J+a^-^=a^-^ + ff{e'). (187) 

Here, the vector Y is given by y — C£'^ + n°- . Considering e ^ 1, we have 
/C^^-* sa "fih"* a-iid if^P ~ '^ib''- With the same discussion, we get 

.%bX"'Y'' = C^Xbti'' - Xbn-'n'' = -Xbu'^n'' + ^(e*) . (188) 

Since the null vector is antiself dual, i.e., eabn^ = ^^a, there is a sign 
difference from the case of the future trapping horizon. Additionally, if we 
assume the dominant energy condition is satisfied, then, we have 

Ika'^fccnll ^ , h^q'^n^^^ebcdW ^ ^ . (189) 
Similar to the future case, the slowly evolving conditions are given: 
(P-i'). \\£x^a\\ and ILxKxl ^ e/L^; 
(P-ii'). ICxO^^^ ^ e/L\ 

With these conditions, one can find that kx is nearly a constant on the past 
trapping horizon. So it can also be expanded as kx — + ff{e). 

At the end of this subsection, some discussions are given as follows: 
(i). In the discussion of the null trapping horizons, for some given foli- 
ation structure (some given X) of the trapping horizon, we can always get 
X^'VfjX" = ±KxX'^ by selecting a special null frame. So, in some sense, we 
can foliate the null trapping horizon arbitrarily. However, in the non-null cases, 
the nonvanished functions C"s carry the information of the foliation, and any 
restriction on them is actually a kind of restriction on the foliations. So, to 
discuss the slowly evolving trapping horizon, we have to select an appropriate 
foliation, i.e., the vector X. 
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(ii) . For some special cases in which the codiniension-2 surface has enough 
symmetry, for example, the spherical symmetry, the problem is greatly simpli- 
fied j50j. In this case, some quantities such as the shear tensor C^j,'^ and the 
5*0(1, 1) connection uja are both vanished. Further, most of the tensors are 
independent of the points on the codimension-2 surface. 

(iii) . Clausiu like equations: For the future slowly evolving trapping horizon, 
substituting the results (|177p . (|178p and (I183|) into eq. ([75)) . and considering 
slowly evolving condition (F-ii'), we get 



\8ttG 



(190) 



Similarly, for the past slowly evolving horizon, we have 



SttG 



Cx.b/ = / e, 



(191) 



These are Clausius like equations. The integrals of the matter flux S/'abi"'^ 



or 3^abn°"n and gravitational radiation \\a 



W|12 



have the form of 



±TCxS. Here S = £//4G and T = \ko\/2tt. Certainly, these kind of Clausius 
like equations hold up to the second order of the e. It should be noted here: On 
the future trapping horizon, we have Cx-s^ — — C J eqO'-^^\ while on the past 
trapping horizon, we have Cx£/ = CJ eq6'(^) . So the sign oi Cx-B^ is determined 
by the function C and the type of the trapping horizon. Assume the null energy 
condition is satisfied, then, the positive temperature requires that SxS > 0. 
Here, we have defined Sx = ±sigu{Ko)Cx , and " + " and " — " correspond to 
the future and past respectively. 



8 Trapping Horizon in FLRW Universe 

In this section, as an example, we study the slowly evolving trapping horizons 
in the FLRW universe. We also discuss the dynamics of these kinds of trapping 
horizons in the formalism with the quasilocal energy (the Misner-Sharp energy) 
at the end of this section. 

8.1 Classification of the Trapping Horizons in FLRW Uni- 
verse 

Firstly, we gives the classification of the trapping horizons in the FLRW universe. 
The metric of the FLRW universe {Ai , g) is 

g = -dt^ + ^—j-^dr^ + a^r^dnl_2 , (192) 

where a — a{t) is scale factor and fc = 0, ±1, while d^l^_2 is the line element of 
an (n — 2)-dimensional sphere. We can decompose the metric into the form as 
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eq.((T]) by introducing two null vectors £ and 



n 



Ld.^^^l(^^dt + ^=^dr^ , (193) 

n^d.^^^-i^-dt-^=L=dr^ . (194) 

So we have hab = ~£an-b — na£h, while qab is just the metric for the sphere part, 
i.e., 

qahdx'^dx^ = a^r^drj^_2. 

Obviously, i and n are both future directed. It's also easy to find £a^° = 
UaTi'^ = 0, £q7i" = —1 and qab^'^ = Qabn"" = 0. Of course, there are some freedom 
to choose these two null vectors, for example, £ \£ and n — > X^^n for some 
positive function A. 

From now on, we will only consider the more interesting case of the four 
dimension. After a simple calculation, the expansions of the sphere along these 
two null directions are given by 



q'^'VaUb ^V2{H-J^-\] . (196) 



.^2 

Here f is defined as f = ar. It's also easy to find 



a 



"9 \ -9 9 



From the expansions in eas. (|195p and (|196p . it's easy to find that: when H < 0, 
we always have 6'*^"'^ < 0. So the trapping horizon is given by 6'(^) = 0, and this 
implies relations 



After substituting above relations, ea. (ll97p becomes 

= ij + 4 , £„0(") = ij - _ A . (199) 
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These mean that the marginal surface is always future. Further, the outer or 
inner of the marginal surface is determined by the sign oi H + 2H^ + k/a^. 
So the marginal surface is outer ii H — k/a^ < ~2{H^ + k/a^) and inner if 
H — k/a^ > —2{H^ + k/a?). When the null energy condition is satisfied, we 
always have H — kjc? < 0. So the marginal surface may be inner or outer. 
For H > 0, we always have 9^^^ > 0. When 6*'"' = 0, we have 



and now eas. (ll97p becomes 



^--^=H, 1 = ^2 + A, (200) 



a2 



/:^6l(") = H + 2i?2 + 4 , /:„6'(") = H-X- (201) 

So the marginal surface is always past. The marginal surface is outer if iJ — 
k/a^ > -2(ij2 + fc/a2) and inner if H ~ k/a^ < ''2{H^ + k/a^). When the nuU 
energy condition is satisfied, we always have H — k/a? <Q. So there are some 
ranges in which the marginal surface is outer or inner if we do not impose some 
additional energy conditions. 

It's easy to find, in the cases with if = 0, we always have 6**^^^ > and 
6l(") < 0, so there are no trapping horizons in this case. This is expectable 
because now the spacetime is actually a flat spacetime. 

The FLRW universe is a typical spherically symmetric spacetime, so, as 
discussed before, the components of the evolution vector X, i.e., A and are 
functions which only depend on the coordinates t and r. So they are constants 
on the marginal surface. It's easy to find, for _ff < 0, on the future trapping 
horizon, the relation CxO'^^^ — gives 

ACtO^'^^ = BCn9^'^^ . (202) 

With the null energy condition, from ea. (l94l) . we always have LiO'^^'^ < 0. Thus, 
for the outer marginal surface, we have Lnd''^^ < 0, then, A and B have same 
signs. While for the inner marginal surface, we have Cnd^^^ > 0, ^ and B have 
opposite signs. In both cases, the relation between A and B is given by 

In the cases > 0, we have past trapping horizon. It's easy to find the 
relation Cxd^n) = on the trapping horizon gives 

ACe9^'''> = BCn9^'''> . (204) 

By using ea. (|94l) . with the nuU energy condition, we always have LnO^^^ < 0. 
So, for the outer marginal surface, we have Ci9^^^^ > 0, and then A and B have 
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opposite signs. While for the inner marginal surface, we get CiO^"'' < 0, so A 
and B have same signs. In both cases, A is given by 

^/ H-k/a' \ 

Conclusively, in FLRW universe, for the future trapping horizon, we can set 
Xa = A{£a — aua) With a given by 

H-k/a^ 

a = -. . (206) 

H + 2m + k/a? 

For the past trapping horizon, X has form Xa = B{aia — no) with the same a 
given in above equation. The classification of "outer" and "inner" are indicated 
by the sign of a. 

The future trapping horizon is outer when a > 0, and inner when a < 0. 
Since in the future cases, we have X°'Xa = 2A^a, so the future outer trapping 
horizon is a spacelike hypersurface, while the future inner trapping horizon is a 
timelike hypersurface. 

Similarly, the past trapping horizon is outer when a < and inner when 
a > 0. In this case, we have X'^Xa — 2B^a, so the past outer trapping horizon 
is timelike, while the past inner trapping horizon is spacelike. 

We have to mention some special cases: For the future trapping horizon, if 
£f9'^^) = 0, then from eq. pl^ . we have a = 0. The trapping horizon reduces to 
a null hypersurface. Similarly, the past trapping horizon also reduces to a null 
hypersurface when a = 0. 

So, for the future and the past trapping horizons, the general form of the 
vector X° is given by — an") and B{a£'^ — n°) respectively. Certainly, 
different A^s and B'a correspond to the different foliations of the trapping hori- 
zons. However, it should be noted here: A and B are both constants on the 
codimension-2 surface. So, for the future trapping horizon, we can view the 
factor A as a factor / provided by some relabeling of the foliation. Similarly, for 
the past trapping horizon, B can be viewed as a relabeling factor by selecting 
some foliation parameter t. By this consideration, the evolution vector X"" has 
the form 

X" = /(r - an") (207) 

in the future case, and 

X"" = f{ar - n") (208) 

in the past case. The null vector still has the form (|193|) and ()194p . In the 
following discussion, we will fix / to be unit. We will find this selection of / 
consists with the requirement |C| — di ^ ■ Here, a plays the role of AB in 

ea. lfTTi]) . 
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8.2 Null Trapping Horizons in FLRW Universe 

For the null trapping horizon, the function a is identically vanished. From the 
expression of a in (|206p , we have 

H -k/a^ =0 

on the trapping horizons. The surface gravity associated with for the null 
future trapping horizon is given by 

= ^ . (209) 

The condition that kx (now X" = £") is a constant on the null future trapping 
horizon gives CiKg — H /2 — (So we can only consider the case with k — 0). 
This means H is a negative constant on the horizon. It's also easy to find 

6l("' = 2V2H , £„6'(^) = 2H^ , de^''^ = , /:„6'(") = -2H^ . (210) 

So there are only null future inner trapping horizon. 

Similarly, for the past case, the surface gravity is given by 



H 

71' 



(211) 



and Cn^n = also gives that H — 0. So is a positive constant on the past 
null trapping horizon. We also have 

6i(^) = 2V2iJ , /If^^") = 2ij2 , CnO^.'^'' = , Ci9^'^'> ^ -2H^ . (212) 

This means there are only null past outer trapping horizon. 

Conclusively, on the null trapping horizons (future and past), the Hubble 
parameter H is always a constant. These kinds of horizons exist only when 
fc = 0. Further, only inner horizon exists in the future case, and only outer 
horizon exists in the past case. In the following discussions, we always set 
/c = 0. It should be noted here: the surface gravity (|209p (or (HH])) is defined 
up to a positive constant coefficient. 



8.3 Slowly Evolving Trapping Horizon in FLRW Universe 

Since the spherical symmetry, it's very simple to study the slowly evolving 
properties of the trapping horizon in the FLRW unverse. In this case, most 
of the scalars on the codimension-2 surface are constants. For example, from 
the definition, the evolving parameter e in the condition (F-i) becomes (we only 
consider the four dimension case, and choose L to be the radius f — l/|i?| for 
yfc^O.) 




(213) 
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Similarly, for the past trapping horizon, the condition (P-i) is given by 



^ = |a| ( 'Sabti' + ) . (214) 




Since / is chosen to be unit, so we have \C\ = \a\. Straightforward calculation 
shows: in both cases, on the trapping horizons, e's are given by 



Remembering that we are considering the trapping horizons which are near the 
null future inner trapping horizon or the null past outer trapping horizon, so 
the functions a's are assumed to be negative. By defining 

■^ = -;^>0, (216) 
then, from the expression of a in eq. (|206p . we have 



s 



2~s 



(217) 



To ensure that a < 0, we have to require s < 2. The evolution parameter e now 
has a simple form 



2 ' 4+s 



It should be emphasized: the result of e in above equation is independent of 
the rescaling of the null frame and the relabeling of the foliation. Under the 
relabeling of the foliation, we have a — > f^a, while the term inside the round 
brackets in eq. (|213|) (or (|214l) ') will accept a factor 1//^. So e is invariant. 

By ea. (l218p . we can express s as a function of e. Considering e <C 1, we have 

So ea. (l217p gives |a| « Actually, from eqs. (|218p and (I217p . we can express 

a as a function of e explicitly. On the other hand, if we do not choose / = 1, 
then we have \C\ = \pa\ w pe^/i. To ensure the relation (|174p . p has to be 
finite and order one. Certainly, the selection of / = 1 satisfies this requirement. 
It's easy to find that 

^ Q Q „ 



Qabe^n" = — + H= . 

So, when s ^ 1, the conditions (F-ii) and (P-ii) are easily satisfied, i.e., Qab^^'n^ ^ 
1/f^ . Obviously, the conditions (F-iii) and (P-iii) are trivially satisfied. Thanks 
to the spherical symmetry and the Weyl flat of the FLRW universe, eqs. (|180p , 



51 



(|18ip and (|189p are also trivially satisfied even without any energy condition 
(In the general case, the dominant energy condition is required). 

For the future inner trapping horizon (H < 0), the slowly evolving conditions 
(F-i') and (F-ii') reduce to \Cxkx\ ^ e/r^ and \CxO^'''^ ^ e/f^. It's easy to 
find 



-H's 1 



2 + .S 
2-s 



(220) 



Therefore, from ea. (l219p . we have 

So, when e <C 1, \Lxd''"'^\ d: ^/r^ is automatically satisfied. After substituting 
Ki = H/\/2 and k„ = we find 



H , 



After a simple calculation, we have 

(2 - sf 



< 0. 







2 - s + + 









Thus, if we require 



H 



< e. 



(221) 



(222) 



(223) 



then we have j/CjfKjfl < e/f^ = eH^ . The condition ([223|) means H is also 
required to be slowly evolving. 

In the past outer case {H > 0), the conditions (P-i') and (P-ii') reduce to 
I'CxKjfl d: f-l'T^ and I^Cx^*^^-*! d e/f"^- A simple calculation shows: on the past 
outer trapping horizon, we have 



2 + 5 



(224) 



So the condition \CxO^^^\ d e/f"^ is also automatically satisfied when e ^ 1. 
Now, it's easy to find the surface gravity kx is given by 



KX 



H 

71V 



2-s 



>0, 



and 







+ + 









(225) 



(226) 



So the condition on |£x'«x| gives same constraint as the one of the future inner 
trapping horizon, i.e. \H\ < eH^- 
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Conclusively, the requirement of the evolving parameter e ^ 1 automatically 
implies that s — —H/H^ is very small. While the slowly evolving condition of 
Kx requires that \H/H^\ is also a small quantity. 

For the slowly evolving past trapping horizon {H > 0), s ~ —H/H^ is very 
small. This is just one of the condition of the slow-roll inflation. Additionally, 
the slowly evolving condition also requires that \H/H^\ is small to ensure that 
the surface gravity kx changes slowly on the trapping horizon. This way, the 
system is near an equilibrium state. This requirement may has some relation 
to the second condition of the slow-roll inflation scenario. Of course, here, we 
have not introduce any scalar field and the corresponding potential. 

For the past horizon, from ea. (|191l) . we have 



^^^CxJ^ ^ J egXhU^'n" . (227) 

Up to second order of e (or the first order of s). Actually, it's easy to find 

Cxs^ = ~V2 (I) (2^) « -^/2 (I) . + ^(.^) (228) 

So we have Lx^i < 0. This just means that £/ decreases along X direction. 
Noted that X° = — n° if s = 0, so X is past pointing. By this consideration, the 
negative Hx^ just means the area of the marginal sphere of the past trapping 
horizon increases along the future direction. The leading order of nx is Hj \/2, 
so we have —{ko/9>ttG)Cx^ = s/2G + ff{s'^). It's also easy to find 

eq^Ta^n'^n^ = ^ . (229) 

Thus, eg. dlQip holds up to the second order of e. It should be emphasized here: 
the surface gravity kx is defined up to a constant factor. Actually, the factor 
•\/2 in Kx comes from the selection of the null vectors (jl93|) and (|194|) . From 
this Clausius relation like equation, we have 

T=^^^fl-f)+^(.4). 



27r 2tt \ 2 

8.4 Quasilocal Energy and Horizon Dynamics 

To compare with the results for the slowly evolving trapping horizon in the 
FLRW universe, in this subsection, we also consider the case where /c = 0. The 
quasilocal energy (|106l) (or the energy (11121) ) inside the sphere with the radius 
f is 

= ^ (1 - V,r~V"f) . (230) 

It's easy to find the effective surface gravity R in ea. (|125p is given by k = l/2f = 
\H\/2 which has the same form of the surface gravity of Schwarzschild black hole. 
While the surface gravity n in ea. (|127p becomes 
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where s is defined in cq. (|216p . So the temperature of the past outer trapping 
horizon is T = |K|/27r = H {1 - s/2) /27r. 

If we omit the factors \/2 in eqs. ()22ip and (|225p . the temperatures corre- 
sponding to K and Kx are coincided up to the second order of e. In the descrip- 
tion with the quasilocal energy, the Clausius relation is given by (corresponding 
to ea. l^n^i ) 

A^aX^^^CxS, (232) 

where S* = / AG and X is the evolution vector of the trapping horizon. Of 
course, above equation is exactly held on the trapping horizon. More detailed 
discussions of the thermodynamics of the FLRW universe in various gravity 
theories can be found in [311 EZl HI [Ml EOl IMl [Ml IMl ES] and related topics 
can be found in a review paper |66j . 



9 Conclusion and Discussion 

In this paper, we have studied the deformation of some spacelike submani- 
fold with an arbitrary codimension. By requiring the projection operator is 
Lie dragging along a given normal vector X, we get ea.(|5(I)). In the case of 
codimension- 1, it reduces to the evolution equation of the extrinsic curvature of 
the spacelike hypersurface ((56)) (and ea.(|55|)). Even in this simple case, there 
are also some interest applications on thermodynamics. For example, in static 
case, the evolution equation (|55l) becomes 

By defining the so called Tolman-Komar mass Mk = J tsNMabU°'u^ inside 
some {n — 2)-closed surface S embedded in the hypersurface, one can give some 
discussion on the thermodynamics associated with the surface S. More details 
can be found in [571 [Ml [M]- In the more interesting case of codimension-2, 
eg. ([SO)) reduces to eg. ([SO)) , i.e., 

-^i (i? - KabcK"''"' ~ KaK') ■ {XeY'') 

-Y^D.D^X, - (X^VeY') . 

This result is frame independent, and it reduces to the well known (cross) focus- 
ing equations after selecting a local null frame. The deformation of the 5*0(1, 1) 
connection is given by 

with Ya = £abX'' if we define an appropriate "surface gravity" kx- It's just 
a Damour-Navier-Stokes like equation if X is self-dual or anti-self-dual. We 
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have investigated the relation between these deformation equations and the 
dynamics of the trapping horizon in two different formahsms: with and without 
introducing some quasilocal energy. In the first formalism, we have proposed a 
generalized energy (|112[) in the higher dimension Einstein gravity theory, i.e.. 



This energy reduces to the Hawking energy in the four dimension. In the case of 
general spherical symmetry, it reduces to the Misner- Sharp energy in the higher 
dimension. We have also studied the deformation of this quasilocal energy, and 
the deformation equation has been given in ea. (|120p . Once the requirement 
CxJ^ = is fulfilled, on the trapping horizon, we always have 



This result is the total variation (evolution) of the energy on the trapping hori- 
zon, and it can be decomposed into two parts as in the four dimension: con- 
tributions from the matter fields (^h) and the contribution from gravitational 
radiation (||cr|p and ||CIP)- When the marginal surface is Einstein, we also study 
the first law like equation of the trapping horizon. It has a similar form as the 
one with the spherical symmetry. However, generally, it's impossible to define 
a surface gravity which is a constant on the marginal surface. Further, it also 
evolves on the trapping horizon even in the spherically symmetric case. This 
means the system is generally nonequilibrium (even far from equilibrium point) 
if we regard the temperature is proportional to the surface gravity. To make the 
problem clear, we have studied some near equilibrium state by considering the 
slowly evolving trapping horizon proposed by Booth in the formalism without 
the quasilocal energy. 

To study the past trapping horizon in the FLRW universe, we generalize 
the definition of the slowly evolving future outer trapping horizon to the past 
trapping horizon. We find, for the slowly evolving past trapping horizon, the 
Clausius like equation is modified to be 



After classifying the trapping horizon in the FLRW universe, as an example, 
we study the slowly evolving trapping horizon in this spacetime. We find, to 
require the past trapping horizon in the FLRW universe to be slowly evolving, 
the Hubble parameter has to satisfy —H /H"^ « ^ 1. Further, \H\/H^ is 
required to be (at most) the order of e. These conditions have close relation to 
the scenario of the slow-roll inflation. We also compare the temperature defined 
in the formalism with the quasilocal energy (T = |/c/27r|) and the temperature 
defined in the formalism without the quasilocal energy ( T = |Kj(-/27r|). We find 
these two temperatures are essentially the same up to the second order of the 
slowly evolution parameter e. 






55 



10 Acknowledgement 



The author would hke to thank Prof. Rong-Gen Cai and Prof. Nobuyoshi Ohta 
for their long term encouragement and kind help. The author would like to 
thank Prof. K.i.Maeda, Prof. S.Mukohyama and Prof. T.Shiromizu for their 
useful discussions on the surface gravity and the generalized energy form. The 
author also thanks Dr.Masato Minamitsuji, Dr.Masato Nozawa and Dr.Chul- 
Moon Yoo for their useful discussions and kind help. The author would like 
to thank Dr. Booth for his valuable suggestions and comments. This work is 
supported by JSPS fellowship No.P09225. 

Appendix A 

In this appendix, we give a detailed derivation of ea.([5(I|. It's easy to find the 
first term in ea.(|49p is given by 

q^^q.'X'^V.KZ^ = 9a'%'^^'=Ve (qJq/VfY,) 

= q,'q,'X-V,qJWfYg + qj q.'^X'^V ,qJ>V fY, + qj q,^ X^V ,V jYg 
= i%'^fYg) DaXf + {qJ^fY,) b,X^ + qJq.^X^V^VfYg 
= i^'^fYs) {DaXf) + [DaY,) 

+ qJq,'X-Y\^,f,^, + qJq^'X^VfV^Yg . (233) 

Here, from second line to third line, we have used ea. (|411) . Again, by using 
eq. (|lT|) or eq. p^ . we have 

-iJlb'^fX^'^eYg - qJWfq,'X^W,Yg 
= - - q/Y,Vf{Di,X') 

-K^V'^Kir - [DaXf) - qJ\/jq,^X'^V,Y, . (234) 

It's also easy to find 

qjY,^f{b,X^) = K^^^'k[^^ + qJqCY^^f^eX, 

+ qJ^f%'y'^aXe. (235) 

So we get 

qJq,'X^VfV,Yg = ~qJVf{Y,D,X^) 

- qJ%'^fX'^eYg ~ qJVfq.'X^WeY, 

- i%'^fYs) {baXf) - qJVfq,^Y^V,X, 

- qJ\/fq,'X^\/,Yg - qJq,^Y^\/f\/,Xg . (236) 
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Then, after substituting above equation into ea. (|233p . we find 

-Ki-)^K^r~K^1''Krj+Kil\ (237) 



where Za = CxYa- Substituting above equations into eg. (1491) . we arrive at 

LxKZ^ = Kif^ + qJq,'X^Y'^^,fgu - qj q.'Y^V fV eXg . (238) 
However, since we have 

qJqb'Y^Vf'^eX, = qjY^S/f {q^V eX,) - qJVfq.'Y^VeX, 

- qJVfq.'^Z, + qJVfq.'^XSVgY, 

= qjY^Vf {Kif^ + DbXg) + q/q.'^VfZ, 

+ q/Vfq,'^{qJ^ + hJ^)X3VgY^ 
= -qJVfY^Kif^ + Y^DaD.X, + K^^ 

+ Kacb (^^VgF") + Kabc (^^VgF") 

= -Ki'-^^Kif^ + Y^DaD.X, + 

- Kacb (YcD'X''^ + Kabc (^^VgF^) , (239) 

and after substituting this resuh into ea. (|238p . we finahy get the deformation 
equation ([SUI- 



Appendix B 

Remembering that eab is defined as 

£ab = naib - (-aUb Or tab = UaVb ~ VaUb , (240) 

it's easy to find 

^a^^cb = hab , ^ab^cd — hadhbc — hachbd ■ (241) 

This also suggests that we can introduce the tensor tab without using the local 
frames. By defining 

j^a ^ ^abj^^ ^ ya ^ ^aby^ ^ (242) 

then, we have 

X" = e^^Xb , Y" = e^^^Yb , XaY" = -XaY" . (243) 

It's also easy to find la — tab^ and Ua = —tab"^- This means: la is self-dual, 
but Ua is anti-self-dual. 
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Appendix C 



From the definition of the covariant derivative D, we have 



— DaX'^eijel + ujabd.X'^ 
^ DaX-'eijel+ujaebdX'^ , 



)Xd 



(244) 



and then 



- {D.DaX'^) eijei + DaX^uj.eMeije''' 

+ DcUJaebdX'^ + DcX-'iOaebdSlje^'^ + UJaL^c£bde'^eX'' ■ 



(245) 



So we get eq. ([6T|) . Here, we have used the fact that Da^tc — 0. 

References 

[1] J. D. Bekenstein, Phys. Rev. D 7, 2333 (1973). 

[2] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [Erratum-ibid. 46, 
206 (1976)]. 

[3] J. M. Bardeen, B. Carter and S. W. Hawking, Commun. Math. Phys. 31, 
161 (1973). 

[4] R. M. Wald, Living Rev. Rel. 4, 6 (2001) | arXiv:gr-qc/9912119| . 

[5] S.W. Hawking, G.F.R. Ellis, The large scale structure of space time, Cam- 
bridge: Cambridge University Press, 1973. 

[6] A. Ashtekar and B . Krishnan, Living Rev. Rel. 7, 10 (2004) 
|arXiv:gr-qc/040704^ . 

[7] T. Jacobson, Phys. Rev. Lett. 75, 1260 (1995) [arXiv:gr-qc/9504004) . 

[8] C. Ehng, R. Guedens, and T. Jacobson, Phys. Rev. Lett. 96, 121301 
(2006). 

[9] S. A. Hayward, Phys. Rev. D 49, 6467 (1994). 
[10] S. A. Hayward, Phys. Rev. D 53, 1938 (1996) |arXiv:gr-qc/9408002| . 
[11] S. A. Hayward, Class. Quant. Grav. 15, 3147 (1998) |arXiv:gr-qc/9710089] . 
[12] R. Penrose, Phys. Rev. Lett. 14, 57 (1965). 



58 



[13] A. Ashtekar, C. Beetle and S. Fairhurst, Class. Quant. Grav. 16, LI (1999) 
[arXiv:gr-qc/9812065 . 

[14] A. Ashtekar, A. Corichi and K. Krasnov, Adv. Theor. Math. Phys. 3, 419 
(2000) arXiv:gr-qc/99050 89. 

[15] A. Ashtekar, C. Beetle and S. Fairhurst, Class. Quant. Grav. 17, 253 (2000) 
|arXiv:gr-qc/990706^ . 

[16] A. Ashtekar, S. Fairh urst and B. Krishnan, Phys. Rev. D 62, 104025 (2000) 
|arXiv:gr- qc/0005083 . 

[17] A. Ashtekar, C. Beetle, O. Dreyer, S. Fairhurst, B. Krishnan, 
J. Lewandowski and J. Wisniewski, Phys. Rev. Lett. 85, 3564 (2000) 
jarXiv:gr-qc/0006006. 

[18] A. Ashtekar and B Krishnan, Phys. Rev. Lett. 89, 261101 (2002) 
|arXiv:gr-qc/020708^ . 

[19] A. Ashtekar and B. Krishnan, Phys. Rev. D 68, 104030 (2003) 
|arXiv:gr-qc/" 0308033 . 

[20] I. Booth and S. Fairhurst, Phys. Rev. Lett. 92, 011102 (2004) 
[arXiv:gr-qc/0307087^ . 

[21] W. Kavanagh and I. Booth, Phys. Rev. D 74, 044027 (2006) 
| arXiv:gr-qc/0603074 . 

[22] L Booth and S. Fairhurst, Phys. Rev. D 75, 084019 (2007) 
|arXiv:gr-qc/0610032l. 

[23] L. B. Szabados, Class. Quant. Grav. 11, 1833 (1994) [arXiv:gr-qc/9402001] . 

[24] K Andersson, M. Ma rs and W. Simon, Phys. Rev. Lett. 95, 111102 (2005) 
|arXiv:gr qc/0506013] . 

[25] L. Andersson, M. Mars and W. Simon, Adv. Theor. Math. Phys. 12, 853 
(2008) larXiv:0704.2889l [gr-qc]. 

[26] C. W. Misner and D. H. Sharp, Phys. Rev. 136, B571 (1964). 

[27] S.W. Hawking, Gravitational radiation in an expanding universe, J. Math. 
Phys. 9, 598 (1968). 

[28] S. A. Hayward, S. Mukohyama and M. C. Ashworth, Phys. Lett. A 256, 
347 (1999) a rXiv:gr-qc/98100 06 . 

[29] S. A. Hayward, Phys. Rev. Lett. 93, 251101 (2004) ^arXiv:gr-qc/040407'i^ . 

[30] S. A. Hayward, Phys. Rev. D 70, 104027 (2004) [arXiv:gr-qc/0408008] . 



59 



[31] H. Bray, S. Hayward, M. Mars and W. Simon, Commun. Math. Phys. 272, 
119 (2007) a rXiv:gr-qc/0603014j . 

[32] R. M. Wald, General Relativity, University of Chicago Press, Chicago, 1984. 

[33] B. Carter, J. Geom. Phys. 8, 53 (1992). 

[34] B. Carter, |arXi7:hep-th/9705172| 

[35] B. Carter, Int. J. Theor. Phys. 40, 2099 (2001) [arXiv:gr-qc/0012036] . 
[36] E. Gourgoulhon,|arXiv:gr-qc/0 703035 | 
[37] T. Damour, Phys. Rev. D 18, 3598 (1978). 

[38] E. Gourgoulhon, Phys. Rev. D 72, 104007 (2005) |arXiv:gr-qc/0508003] . 

[39] E. Gourgoulhon and J. L. JaramiUo, Phys. Rev. D 74, 087502 (2006) 
|arXiv:gr-qc/0607050 . 

[40] E. Gourgoulhon and J. L. Jaramillo, New Astron. Rev. 51, 791 (2008) 
[arXiv:0803.2944 [astro-ph]]. 

[41] J. D. Brown and J. W. York, Phys. Rev. D 47, 1407 (1993) 
|arXiv:gr-qc/920901^ . 

[42] K.S. Throne, R.H.Price, and D.A.MacDonald, Black Holes: The Membrane 
Paradigm, Yale University Press, New Haven, 1986. 

[43] M. Parikh and F . Wilczek, Phys. Rev. D 58, 064011 (1998) 
|arXiv :gr-qc/9712077l. 

[44] C. Eling, I. Fouxon and Y. Oz, Phys. Lett. B 680, 496 (2009) 
^Xiv:0905.3638i [hep-th]]. 

[45] C. Eling and Y. Oz, JHEP 1002, 069 (2010) |arXiv:0906.4999l [hep-th]]. 

[46] I. Booth, M. P. Heller and M. Spalinski, Phys. Rev. D 80, 126013 (2009) 
|arXiv:0910.0748l [hep-th]]. 

[47] R. Mizuno, S. Ohashi and T. Shiromizu, Phys. Rev. D 81, 044030 (2010) 
^arXiv:0911.5560 [gr-qc]]. 

[48] H. Kodama, Prog. Theor. Phys. 63, 1217 (1980). 

[49] A. B. Nielsen and J. H. Yoon, Class. Quant. Grav. 25, 085010 (2008) 
[a rXiv:071 1. 1445 [gr-qc]]. 

[50] I. Booth and J. Martin. a rXiv: 1007. 164 2 [gr-qc]. 

[51] H. Maeda and M. Nozawa, Phys. Rev. D 77, 064031 (2008) 
[arXiv:0709.1 199 [hep-th]]. 



60 



[52] M. Nozawa and H. Maeda, Class. Quant. Grav. 25, 055009 (2008) 
^arXiv:0710.2709 [gr-qc]]. 

[53] A. L. Besse, Einstein Manifolds, Springer- Ver lag, 1987. 

[54] S. A. Hayward, Phys. Rev. D 74, 104013 (2006) [arXiv:gr-qc/0609008] . 

[55] S. A. Hayward, Class. Quant. Grav. 24, 923 (2007) [ arXiv:gr-qc/0611027] . 

[56] A. Frolov and L. Kofman, JCAP 0305, 009 (2003) 

|arXiv:hep-tli/021232 7 . 

[57] U. H. Danielsson, Phys. Rev. D 71, 023516 (2005) |arXiv:hep-th/0411172] . 

[58] R. Bousso, Phys. Rev. D 71, 064024 (2005) [arXiv:hep-th/0412197] . 

[59] G. Calcagni, JHEP 0509, 060 (2005) |a rXiv:hep-th/0507125| . 

[60] R. G. Cai and S. P. Kim, JHEP 0502, 050 (2005) [arXiv:hep-th/0501055] . 

[61] R. G. Cai, L. M. Cao and N. Ohta, Phys. Rev. D 81, 084012 (2010) 
[arXiv:1002.1136 [hep-th]]. 

[62] R. G. Cai, L. M. Cao and Y. P. Hu, JHEP0808: 090,2008 [arXiv: 
0807.1232[he p-th]]; M . Akbar and R. G. Cai, Phys. Rev. D 75, 084003 
(2007) |arXiv: hep-th/ 0609128 ; R. G. Cai and L. M. Cao, Phys. Rev. D 75, 
064008 (2007)"]aFXw:gr-qc/061107f ; Y. Gong and A. Wang, Phys. Rev. 
Lett. 99, 211301 (2007) arXiv:0704.0793 [hep-th]]; S. F. Wu, B. Wang and 
G. H. Yang, Nucl. Phys. B 799, 330 (2008) ^arXi^07lT7l209l [hep-th]] . 

[63] R. G. Cai, P rog. Theor. Phys. Suppl. 172: 100-109, 2008 
[arXiv:0712.2142;hep-th]]; R. G. Cai and L. M. Cao, Nucl. Phys. B 785, 
135 (2007) |arXiv:hep-th/06 12144 ; A. Sheykhi, B. Wang and R. G. Cai, 
Nucl. Phys. B 779, 1 (2007) narXiv:hep-th/0701198]; A. Shey khi, B. Wang 
and R. G. Cai, Phys. Rev. D 76, 023515 (2007) 'arXiv:hep-th/ 0701261] ; 
X. H. Ge, Phys. Lett. B 651, 49 (2007) |arXiv:hep-th/0703253 . 

[64] T. Padmanabhan, Class. Quant. Grav. 19, 5387 (2002) 
[arXiv:gr-qc /0204019 ] ; D. Kothawala, S. Sarkar and T. Padmanab- 
han, Phys. Lett. B 652, 338 (2007) arXiv:gr-qc/0701002 ; A. Paranjape, 
S. Sarkar and T. Pa dmanabhan, Phys. Rev. D 74, 104015 (2006) 
arXiv:hep-th/060 7240| ; R. G. Cai, L. M. Cao, Y. P. Hu and S. P. Kim, 
Phys. Rev. D 787T24012 (2008) arXiv:0810.2610 [hep-th]]; M. Akbar, 
Chin. Phys. Lett. 24, 1158 (2007) arXiv:hep-th/070202 9 ; M. Akbar and 
A. A. Siddiqui, Phys. Lett. B 656, 21^(2007); M. Akb ar and R. G. Cai, 
Phys. Lett. B 648, 243 (2007) jarXiv:gr-qc/0612089| ; R. G. Cai and 
N. Ohta, arXiv:0910.2307 [hep-th]. 

[65] K. i. Maeda and M. Nozawa, Phys. Rev. D 81, 044017 (2010) 
|arXiv:0912;28TTl [hep-th]]: K. i. Maeda and M. Nozawa, Phys. Rev. D 81, 
124038 (2010) [arXiv: 1003.28491 [gr-ac]]. 



61 



[66] T. Padmanabhan. larXiv:091 1.5004 ' [gr-qc]. 

[67] T. Padmanabhan, Phys. Rev. D 81, 124040 (2010) [arXiv: 1003.56651 [gr- 
qc]]. 

[68] E. P. Verlinde. larXiv: 1001. 0785] [hep-th]. 

[69] G. Abreu and M. Visser, Phys. Rev. Lett. 105, 041302 (2010) 
|arXiv:1005. 11321 [gr-qc]]. 



62 



